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ABSTRACT 

A 

Far-field  Rayleigh  and  Love  waves  observed  on  a long-period 

seismic  network  were  studied  to  characterize  presumed  underground  nuclear 
explosions.  These  occurred  at  the  United  States  Nevada  Test  Site  (NTS),  the 
Russian  eastern  Kazakh  (EKZ)  region,  and  at  the  site  of  presumed  peaceful 
nuclear  explosions  (PNE)  in  the  north  Caspian  Sea  region.  The  observation 
stations  available  for  this  study  included  Alaskan  Long  Period  Array  (ALPA), 
Large  Aperture  Seismic  Array  (LASA),  Norwegian  Seismic  Array  (NORSAR), 
the  Very  Long  Period  Experiment  (VLPE)  stations,  the  Seismic  Research  Ob- 
servatories (SRO),  and  the  Special  Data  Collection  System  (SDCS)  stations. 

The  purpose  of  this  study  was  to  determine  the  seismic  source  parameters 
from  long-period  surface  wave  data  and  to  investigate  the  influence  of  geo- 
logical and  near- source  propagation  effects  on  the  frequency  content  of  the 
surface  wave  signals. 

The  surface  wave  group  velocities  were  measured  for  travel 
paths  from  NTS,  EKZ,  and  PNE  to  all  available  observation  stations.  Results 
were  examined  using  all  selected  NTS,  EKZ,  and  PNE  events.  The  collective 
behavior  of  the  surface  wave  group  velocities  for  each  travel  path  were  then 
obtained. 

Amplitude  spectra  of  the  NTS,  EKZ,  and  PNE  events  under 
study  were  examined  in  the  period  range  between  10  and  50  seconds.  The  am- 
plitude spectral  ratios  of  the  Love  wave  to  Rayleigh  wave  (LQ/LR)  were  cal- 
culated for  each  event  at  all  available  observation  stations.  These  ratios 
were  then  compared  among  the  NTS,  EKZ,  and  PNE  events  for  any  indications 
of  geological  and/or  source-oriented  dependence,  such  as  the  correlation  be- 
tween the  LQ/LR  ratio  and  the  event  m , and  any  apparent  differences  in  the 
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LQ/LR  ratios  among  the  NTS,  EKZ,  and  PNE  events,  or  between  the  source 
types,  i.  e.  , explosion  versus  earthquake.  Also,  the  relative  excitations  of 
the  Love  waves  for  a suite  of  NTS  events  were  examined. 


The  average  surface  wave  attenuation  coefficients  along  the 
travel  paths  between  the  NTS  and  the  observation  stations  and  those  along  the 
travel  paths  between  the  EKZ  and  the  observation  stations  were  estimated  by 
the  isotropic- source  method.  All  available  events  were  used.  The  travel 
paths  between  the  NTS  and  the  observation  stations  were  divided  into  two 
groups;  purely  oceanic  path  and  mostly  continental  path.  In  addition,  when- 
ever applicable,  the  two-sta^^ion  method  was  employed  to  calculate  the  surface 
wave  attenuation  coefficients  for  the  travel  path  between  two  observation  sta- 
tions. These  results  were  compared  with  the  existing  results  on  the  surface 
wave  attenuation. 

Finally,  the  seismic  source  parameters  were  estimated  for  all 
selected  NTS,  EKZ,  and  PNE  events  based  on  the  observed  long-period  sur- 
face wave  spectra.  The  estimated  seismic  source  parameters  included  (1) 
relative  strength  of  the  explosive  source  and  double-couple  source  which 
modeled  the  possible  tectonic  strain  release  associated  with  the  explosion, 

(2)  focal  depth,  and  (3)  the  source  mechanism  of  the  possible  tectonic  strain 
release,  i.  e.  , the  dip,  slip,  and  strike  angle.  Layered  half- space  earth 
models  were  obtained  for  the  NTS,  the  EKZ,  and  the  north  Caspian  Sea  re- 
gion. This,  together  with  the  theoretical  source  model  consisting  of  point 
explosive  source  and  point  double-couple  source,  was  used  to  generate  theo- 
retical surface  wave  spectra.  The  results  from  the  source  parameter  esti- 


mations were  then  used  in  a comparative  study  of  NTS,  EKZ,  and  PNE  events 
for  possible  application  to  discrimination  of  earthquakes  from  explosions. 
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SECTION  I 
INTRODUCTION 


This  report  presents  the  results  from  a continuing  effort  at 
Texas  Instruments  Incorporated  to  obtain  the  seismic  source  parameter  esti- 
mates and  to  investigate  the  travel  path  effects  using  long- period  teleseismic 
surface  wave  data.  Previous  work  (Tsai,  1972a,  1972b;  Tsai  and  Shen,  1972; 
Turnbull, et  al,  , 1973,  1974,  1975;  Sun,  1976)  has  been  mainly  concerned  with 
methodology  (such  as  the  development  of  multi-period  multi- site  amplitude 
spectral  fitting  process),  study  of  the  theoretical  surface  wave  spectra  and 
radiation  patterns,  and  verification  and  application  of  the  methods  to  synthetic 
events  and  earthquakes.  All  these  studies  have  been  aimed  toward  the  better 
understanding  and  utilization  of  the  long- period  teleseismic  surface  waves  for 
source  characterization  and  possibly  for  source  discrimination.  Previously 
developed  methods  of  narrowband  filtering  and  amplitude  speciral  fitting  have 
been  proved  to  be  efficient  and  quite  satisfactory  for  the  extraction  of  the  de- 
multipathed  seismic  surface  wave  spectra  and  the  determination  of  seismic 
source  parameters  from  the  far- field  seismic  surface  waves  of  earthquakes. 

In  this  report,  we  apply  the  same  techniques,  making  necessary  modifications, 
to  underground  nuclear  explosions  for  similar  purposes,  namely  the  deter- 
mination of  the  seismic  source  parameters  from  long-period  seismic  surface 
wave  data  and  investigation  of  the  influence  of  the  geological  and  near- source 
propagation  effects  on  the  frequency  content  of  these  surface  waves. 

Since  Tsai's  method  was  established,  many  seismic  events 
have  been  analyzed  by  several  authors  using  this  method.  An  error  analysis 
was  attempted  during  their  application  of  this  technique.  However,  their  en- 
deavor was  defeated  by  many  unsolvable  difficulties.  First,  there  was  no 
standard  solution  which  they  could  compare  their  result  with,  so  no  definite 
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error  could  be  estimated.  Secondly,  errors  for  a system  depend  upon  basic 
assumptions,  the  procedures  used,  the  technique  of  computation,  and  ana- 
lyst's judgments.  There  are  some  high  degree  of  uncertainties  such  as  the 
earth  model  assumptions,  travel  path  attenuations,  data  processing  etc. 

Most  of  these  uncertainties  are  complexly  interrelated.  This  causes  the  dif- 
ferent type  of  errors  to  be  interlocked  in  a way  difficult  to  analyze.  There- 
fore, it  is  not  now  known  by  us  how  to  calculate  individual  error  of  each  step 
of  processing  and  to  accurately  estimate  confidence  limits  of  our  resultant 
source  mechanism  solutions. 

Even  though  the  extent  of  the  errors  could  not  be  determined, 
the  authors  still  thought  this  method  reliable  after  comparing  their  results 
with  those  obtained  by  other  means. 

Battis  and  Hill  (1977)  used  this  method  to  analyze  the  Pocatello 
Valley  earthquake.  The  results  were  similar  to  the  composite  fault  based  on 
a microearthquake  study  of  the  source  region  by  Smith  and  Shar  (1974),  and 
also  agreed  closely  with  that  obtained  by  Arabasz,  et  al.  (1975)  and  Rodgers 
(1977,  from  T.  C.  Bache's  report).  Battis  and  Hill  (1977)  also  calculated  the 
fault  plane  by  this  method  for  the  main  shock  and  the  foreshock  of  Yellowstone 
earthquake.  Their  solutions  are  nearly  identical  to  a composite  fault  plane 
constructed  for  the  same  earthquake  and  its  aftershocks  by  the  U.  S.  G.  S. 
based  on  data  from  a local  array  (C.  Weaver,  1976,  personal  communication). 
Turnbull,  Sun,  and  Battis  (1975)  analyzed  a Bear  Valley  event.  The  mech- 
anism of  that  event  was  reasonably  close  to  the  solution  estimated  by  other 
means.  After  using  this  method  to  analyze  seven  Eurasian  events,  Turnbull, 
Sun,  and  Black  (1974)  found  that  depth  estimates  can  be  obtained  which  closely 
agree  with  those  determined  from  bodywave  depth  phases  (see  Table  I-l). 

Previous  studies  on  the  determination  of  source  parameters  of 
underground  nuclear  explosions  using  far- field  surface  waves  have  been  re- 
ported by  several  authors  for  various  purposes  (Toksoz,  et  al.  , 1964,  1965; 
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Lambert,  et  al.  , 1972;  Mitchell,  1975).  However,  there  are  several  common 
features  in  all  of  these  studies:  (1  ) they  analyzed  a small  number  of  events 
from  the  United  States  Nevada  Test  Site,  typically  two  or  three  explosions, 

(2)  they  used  radiation  pattern  fitting  requiring  many  well  distributed  observa- 
tion stations,  (3)  they  assumed  that  explosions  took  place  at  the  surface  and 
that  the  source  mechanism  of  the  explosion-associated  tectonic  strain  release 
was  identical  to  that  of  the  local  earthquakes  in  the  region,  and  (4)  hence,  to 
obtain  the  best  fit  between  the  observed  and  theoretical  radiation  patterns, 
they  only  estimated  two  parameters  - the  strike  angle,  adjusting  it  around 
the  vicinity  of  the  strike  direction  of  local  earthquakes;  and  the  relative 
strength  of  the  explosion  and  explosion-associated  tectonic  strain  release  i.  e.  , 
the  so-called  F value.  In  this  report,  we  analyze  a total  of  twenty-seven  under- 
ground nuclear  explosions  from  three  different  regions,  the  United  States  Ne- 
vada Test  Site  (NTS),  the  Russian  eastern  Kazakh  (EKZ),  and  the  north  Caspian 
Sea  area.  Also,  we  estimate  all  six  source  parameters  of  a given  underground 
explosion,  namely  the  F value,  focal  depth,  dip  angle,  slip  angle,  strike  angle, 
and  seismic  moment  of  the  explosion-associated  tectonic  strain  release. 

The  following  sections  present  the  results  of  research  on  long- 
period  teleseismic  surface  waves  of  twenty- seven  NTS,  EKZ,  and  PNE  events. 
In  Section  II,  we  discuss  the  data  base  used  in  this  study,  such  as  data  avail- 
ability, data  quality,  and  the  data  reduction  procedure.  The  results  from  the 
examination  of  the  observed  surface  waves  are  presented  in  Section  III,  where 
the  surface  wave  group  velocities  for  various  travel  paths,  the  observed  am- 
plitude spectra,  and  spectral  ratios  are  given  and  discussed.  Section  IV 
briefly  describes  the  methods  which  are  used  for  the  estimation  of  the  surface 
wave  travel  path  attenuations  and  shows  the  attenuation  coefficients  along  var- 
ious travel  paths  obtained  from  this  study.  The  source  parameter  estimations 
of  the  selected  explosions  are  detailed  in  Section  V.  Finally,  in  Section  VI,  a 
summary  of  the  results  given  in  the  previous  sections  is  presented. 


SECTION  II 


DATA  BASE 

A.  SELECTION  OF  THE  EVENTS 

In  the  selection  of  the  events  to  be  analyzed  in  this  study,  we 
tried  to  include  as  many  recent  NTS,  EKZ,  and  PNE  events  as  possible.  The 
NTS  and  the  EKZ  events  are  underground  nuclear  explosions  from  the  United 
States  Nevada  Test  Site  and  from  the  Russian  eastern  Kazakh,  respectively, 
while  the  PNE  events  are  peaceful  nuclear  explosions  in  the  Russian  Eurasian 
continent  some  distance  away  from  the  normal  Russian  test  sites.  A list  of 
the  events  was  obtained  for  the  NTS  and  EKZ  events  by  means  of  a chronolog- 
ical search  (starting  from  1976  and  progressing  backward  into  1975,  1974, 
and  so  on)  through  the  NORSAR  (Norwegian  Seismic  Array)  bulletin,  and  the 
PDE  (Preliminary  Determination  of  Epicenters)  bulletin.  For  the  PNE  events, 
the  list  was  obtained  mainly  from  the  list  of  the  Eurasian  PNE  events  given  by 
Nordyke  (1974).  The  preliminary  list  of  events  consisted  of  twenty-four  NTS 
events  (ten  in  1976,  eleven  in  1975,  and  three  in  1974),  twenty-five  EKZ  events 
(five  in  1976,  seven  in  1975,  seven  in  1974,  and  six  in  1973)  and  ten  PNE 
events  (one  in  1976,  four  in  1973,  three  in  1972,  and  two  in  1971).  Based  on 
this  preliminary  list,  the  availability  of  the  long- period  teleseismic  surface 
wave  data  with  good  quality  was  examined  and  the  final  list  of  the  selected 
NTS,  EKZ,  and  PNE  events,  as  given  in  Table  II- 1,  was  obtained. 

B.  DATA  AVAILABILITY  AND  QUALITY 

In  the  work  statement  of  this  task,  we  were  requested  to  use 


the  surface  wave  data  from  the  observation  stations  which  included  the  Very 
Long  Period  Experiment  (VLPE)  Stations,  the  Seismic  Research  Observatories 


(SRO),  Special  Data  Collection  System  (SDCS)  stations,  Large  Aperture  Seis- 
mic Array  (LASA),  Norwegian  Seismic  Array  (NORSAR),  Iranian  Long  Period 
Array  (ILPA),  and  Korean  Seismic  Research  Station  (KSRS).  The  intention  of 
including  as  many  existing  long-period  recording  stations  as  possible  is  under- 
standable - it  provides  better  azimuthal  coverage  and  better  data  availability. 
Unfortunately,  the  ILPA  data  were  not  available,  since  ILPA  started  opera- 
tion on  1 May  1976  and  almost  all  of  the  events  in  our  preliminary  list  occur- 
red before  that  time;  while  the  KSRS  data  were  not  available  because  no  soft- 
ware was  available  for  editing  the  KSRS  long-period  field  tapes.  To  make  up 
for  the  loss  of  data  from  the  ILPA  and  the  KSRS,  we  included  the  Alaskan 
Long  Period  Array  (ALPA), 

After  checking  the  availability  of  the  surface  wave  recordings 
and  examining  the  results  of  the  processed  surface  wave  signals  from  all 
available  recordings,  the  final  list  of  the  selected  NTS,  EKZ,  and  PNE  events 
was  obtained  as  shown  in  Table  II- 1.  Table  II- 1 gives  the  event  name,  date, 
origin  time,  location,  and  PDE  m^,  together  with  the  availability  of  good 
quality  surface  wave  recordings  indicated  by  the  letter  'B',  'R',  or  'L'  under 
the  'Station  I.  D.  ' heading.  Stations  are  identified  by  the  station  number  (1,2, 

. , , ) or  letter  (A,  L,  andN).  'B'  means  that  both  Rayleigh  and  Love  wave  re- 
cordings are  useable,  'R'  means  that  only  the  Rayleigh  wave  recording  is  use- 
able,  and  'L'  means  that  only  the  Love  wave  recording  is  useable.  The  de- 
scriptions of  the  available  observation  stations  are  given  in  Table  U-2,  in- 
cluding station  name,  location,  and  identification  number  or  letter.  It  can  be 
seen  that  the  final  list  of  the  selected  events  consists  of  twenty- seven  events 
(sixteen  NTS  events,  eight  EKZ  events,  and  three  PNE  events)  instead  of 
fifty-nine  events  in  the  preliminary  list.  The  loss  of  the  events  from  the  pre- 
liminary list  was  mainly  caused  by  (1)  bad  field  tapes  which  could  not  be 
edited  for  the  desired  events,  (2)  bad  edited  traces  which  contained  no  ap- 
parent desired  surface  wave  signals  due  to  the  noise  or  contained  uncorrect- 
able  spikes  in  the  signal  gate,  and  (3)  mixed  events.  In  addition,  we 
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TABLE  II-2 


DESCRIPTION  OF  OBSERVATION  STATIONS 


Station  I.  D. 

Station  Location 

Name 

City  and  Country 

Lat. 

°N 

Lon. 

1 

CTA 

Charters  Towers,  Australia 

-20.  I 

146.  3 

2 

CHG 

Chiang  Mai,  Thailand 

18.  8 

99.  0 

tfi 

4 

TLO 

Toledo,  Spain 

39.  0 

-4.  0 

c 

0 

5 

EIL 

Eilat,  Israel 

29.  6 

35.  0 

OJ 

cn 

6 

KON 

Kongsberg,  Norway 

59.7 

9.  6 

w 

Oh 

7 

OGD 

Ogdensburg,  New  Jersey 

41.1 

-74.  6 

> 

8 

KIP 

Kipapa,  Hawaii 

21.4 

-158. 0 

9 

ALQ 

Albuquerque,  New  Mexico 

34.  9 

-106.  5 

10 

ZLP 

La  Paz,  Bolivia 

-16.  5 

-68.  1 

11 

MAT 

Matsushiro,  Japan 

36.  5 

138.  2 

ANMO 

Albuquerque,  New  Mexico 

34.  9 

-106. 5 

0 

6 

GUMO 

Guam,  Marianas  Island 

13.  6 

144.  9 

CJ 

trt 

7 

MAIO 

Mashhad,  Iran 

36.  3 

59.  5 

12 

TATO 

Taipei,  Taiwan 

25.  0 

121. 5 

>s 

^A 

ALPA 

Alaska 

65.  2 

-147.  7 

nj 

in 

L 

LASA 

Montana 

46.  7 

-106.  2 

N 

NORSAR 

Norway 

60.  8 

10.  9 

03 

c 

1 

FN-WV 

Franklin,  N.  Virginia 

38.  6 

-79.  5 

0 

2 

HN-ME 

Houlton,  Maine 

46.  2 

o 

00 

sO 

1 

C/) 

3 

RK-ON 

Red  Lake,  Ontario 

50.  8 

-93.  7 

V) 

U 

Q 

4 

WH2YK 

White  Horse,  Yukon 

60.  7 

-135.  0 

cn 

5 

CPSO 

McMinnville,  Tennessee 

35.  6 

85.6 

discarded  those  events  which  had  less  than  three  good  edited  surface  wave 
traces  at  two  observation  stations.  This  last  condition  for  rejecting  an  event 
was  based  on  the  previous  result  that  the  source  parameters  of  a given  seis- 
mic event  could  be  satisfactorily  estimated  with  good  quality  surface  wave  re- 
cordings at  two  observation  stations  in  two  consecutive  quadrants  (Turnbull, 
et  al.  , 1973).  From  the  data  availability  indicated  in  Table  II- 1,  it  is  clear 
that  the  NTS  events  have  fairly  good  station  coverage  while  for  the  EKZ  and 
the  PNE  .ivents  we  mostly  have  to  rely  on  the  array  stations. 

The  selected  NTS  events  took  place  at  five  different  locations; 
the  EKZ  events  at  three  different  locations.  Table  II- 3 presents  the  station 
azimuthal  angles  with  respect  to  those  event  locations  and  the  epicenter  dis- 
tances from  those  event  locations  to  each  observation  station.  In  general, 
the  station  coverage  for  the  selected  events  is  thought  to  be  satisfactory  in 
terms  of  station  azimuthal  distribution  and  epicenter  distances  for  the  pur- 
poses of  this  study. 

C.  DATA  HANDLING 

The  observed  surface  wave  data  needed  in  the  analysis  of  this 
study  are  the  amplitude  spectra  from  ten  to  fifty  seconds  for  the  following 
three  purposes: 

• Examination  of  the  surface  wave  spectral  contents  and  spectral 

ratios, 

• Estimation  of  the  surface  wave  travel  path  attenuation,  and 

• Estimation  of  the  seismic  source  parameters. 

In  order  to  obtain  the  desired  surface  wave  amplitude  spectra, 
the  field  tapes  containing  the  surface  wave  recordings  for  the  selected  NTS, 
EKZ,  and  PNE  events  at  all  available  observation  stations  were  processed  as 
follows : 
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The  desired  surface  wave  recordings  on  the  field  tape  are  edit- 
ed to  yield  the  edited  event  signal  which  usually  consists  of 
3«N  signal  traces,  where  N is  the  number  of  recording  sites 
in  a station,  i,  e,  , N = 1 for  a single-site  station.  In  general,  the 
3‘  N signal  traces  in  the  edited  event  signal  are  the  vertical 
and  two  horizontal  (oriented  in  N-S  and  E-W  direction)  record- 
ings, The  vertical  Rayleigh  wave  signal  is  obtained  from  the 
vertical  recording  while  the  radial  Rayleigh  wave  and  trans- 
verse Love  wave  signal  are  obtained  by  properly  rotating  the 
two  horizontal  recordings.  For  the  array  stations,  in  addition 
to  the  rotation  of  horizontal  recordings,  the  edited  event  signal 
is  beamformed. 

The  edited  event  signal  with  proper  rotation  of  its  traces  (and 
beamforming  for  array  stations)  is  reformatted  to  form  the 
event  files  on  7-track  tape  to  drive  the  processing  programs 
in  the  PDP-15/50  computer. 

For  each  event  signal  in  the  event  files,  each  component  trace 
(Rayleigh  or  Love  wave)  is  analyzed  for  evidence  of  multipath- 
ing  by  applying  a series  of  narrowband  filters  (as  discussed  by 
Turnbull,  et  al.  , 1974;  Sun  and  Shaub,  1977).  At  each  period, 
group  velocities  are  determined  for  each  multipath,  and  then 
compared  to  the  standard  group  velocity  curves  for  the  type 
of  travel  path  encountered  (see  example  in  Sun,  1976).  Esti- 
mating the  correct  multipath  from  this  comparison,  the  am- 
plitude spectra  are  obtained.  These  raw  observed  amplitude 
spectra  are  then  corrected  for  the  station  instrument  response 
to  yield  the  desired  observed  amplitude  spectra  for  the  estima- 
tion of  the  surface  wave  travel  path  attenuation.  These  raw 
observed  amplitude  spectra  with  appropriate  instrument 


response  correction  are  further  corrected  for  the  geometrical 
spreading  and  the  travel  path  attenuation  to  yield  the  corrected 
observed  amplitude  spectra  for  the  examination  of  the  surface 
wave  spectra  and  the  estimation  of  the  seismic  source  param- 
eters. 

In  the  above  described  surface  wave  data  reduction  procedure, 
the  first  two  steps  and  the  spectral  correction  in  the  third  step  are  carried 
out  by  batch  processing  on  the  IBM  360/44  computer.  The  third  step  involv- 
ing the  narrowband  filtering  is  done  via  the  Interactive  Seismic  Processing 
System  (ISPS)  developed  by  Texas  Instruments  using  PDP- 15/50  interactive 
graphics. 
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SECTION  III 


EXAMINATION  OF  THE  OBSERVED  SURFACE  WAVES 


Before  applying  the  amplitude  spectral  fitting  technique  to 
selected  NTS,  EKZ,  and  PNE  events,  the  observed  surface  waves  from  these 
events  will  be  first  examined  in  terms  of  (1)  observed  fundamental  mode 
group  velocities  and  spectra,  (2)  amplitude  spectra  and  spectral  ratios  be- 
tween periods  of  10  to  50  seconds,  and  (3)  average  surface  wave  attenuation 
coefficients  estimated  from  observed  amplitude  spectra  along  various  travel 
paths.  The  reasons  for  this  examination,  aside  from  fulfilling  the  work  state 
ment  of  this  task,  are  twofold;  (1 ) the  results  might  be  used  to  guide  future 
processing  and  analysis  of  the  surface  waves  of  other  NTS,  EKZ,  and  PNE 
events,  and  (2)  the  results  might  yield  a better  geophysical  understanding  of 
the  excitation  and  propagation  of  seismic  surface  waves  from  the  NTS,  EKZ, 
and  PNE  events.  The  results  of  the  group  velocity  and  amplitude  spectral 
analyses  are  presented  and  discussed  in  this  section.  The  determination  of 
the  average  surface  wave  travel  path  attenuation  will  be  given  in  Section  IV. 

A,  TRAVEL  PATH  CROUP  VELOCITIES 

In  narrowband  filtering  to  demultipath  surface  wave  spectra  of 
selected  NTS,  EKZ,  and  PNE  events,  group  velocities  were  determined  for 
the  fundamental  mode  great  circle  path  arrival.  This  was  done  for  all  ob- 
servation stations  and  was  repeated  for  all  selected  NTS,  EKZ,  and  PNE 
events.  Therefore,  for  each  event- station  path  (e.  g.  , from  the  NTS  to 
ALPA),  N group  velocity  measurements  were  determined  from  N available 
events  from  the  same  vicinity.  Theoretically,  these  N sets  of  group  velocit 
ies  for  the  same  event- station  path  should  be  identical  since  the  surface  wave 


group  velocities  for  a given  mode  depend  only  on  the  medium  through  which 
the  surface  wave  propagates,  regardless  of  different  seismic  sources.  How- 
ever, in  the  real  world,  there  may  exist  variations  among  these  N sets  of 
group  velocities  caused  by  some  natural  (known  or  unknown)  or  man-made 
factors.  By  examining  these  N sets  of  group  velocities,  some  collective 
behavior  of  the  group  velocities  along  a given  event- station  path  can  be  obtain- 
ed and  might  be  used  for  the  calibration  of  this  particular  event- station  path 
for  future  use. 

Figures  III-l  and  III- 2 show  the  geographical  locations  of  the 
available  observation  stations  in  relation  to  the  NTS  and  the  EKZ  events,  to- 
gether with  the  event- station  travel  paths.  Except  for  the  travel  paths  from 
NTS  to  the  observation  stations  in  the  Pacific  Ocean  which  are  entirely  ocean- 
ic, the  travel  paths  are  either  entirely  continental  or  mixed.  Figures  III- 3 
and  III-4  present  the  Rayleigh  and  Love  wave  group  velocities,  respectively, 
for  the  travel  paths  from  NTS  to  the  observation  stations,  using  all  selected 
NTS  events.  Similar  results  are  shown  in  Figures  III- 5,  IU-6,  III-7,  and 
III-8  for  the  travel  paths  from  EKZ  and  PNE  to  the  observation  stations.  In 
these  figures,  N is  the  number  of  events  used  in  forming  the  average  group 
velocities  for  the  indicated  travel  path,  the  vertical  bar  at  each  given  period 
indicates  the  range  of  the  variation  in  the  group  velocities  obtained  from  these 
N events,  and  the  solid  dot  is  the  mid-point  of  this  range.  Two  solid  curves, 
which  represent  the  theoretical  group  velocity  curves  for  the  typical  continen- 
tal and  oceanic  paths,  are  superimposed  in  these  figures  for  reference  pur- 
poses. 

Referring  to  Figure  III- 3,  which  gives  the  Rayleigh  wave  group 
velocities  for  the  travel  paths  from  NTS  to  the  observation  stations,  the  fol- 
lowing remarks  about  the  behavior  of  these  group  velocities  can  be  made: 

• If  we  take  the  solid  dots  as  the  average  group  velocities  for 

the  indicated  travel  path  and  connect  them  with  straight  lines. 
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a fairly  smooth  group  velocity  curve  can  be  formed,  except  for 
the  travel  path  from  NTS  to  LASA.  Thus  no  sudden  changes  in 
the  velocity  profiles  of  the  earth  along  these  travel  paths  are 
encountered. 

• The  level  of  the  average  group  velocity  curves  can  be  well  ex- 
plained as  typical  continental  or  oceanic  paths,  except  for  the 
travel  paths  to  the  observation  stations  which  are  to  the  east 
of  the  NTS  on  the  northern  American  continent. 

• For  the  travel  paths  to  the  observation  stations  at  LASA,  OGD, 
ALQ,  FN-WV,  HN-ME,  RK-ON,  and  CPSO  which  are  on  the 
northern  American  continent  and  east  of  NTS,  the  levels  of  the 
average  group  velocity  curves  are  lower  than  that  for  the  typi- 
cal continental  path,  especially  at  the  longer  periods  and  for 
the  travel  paths  to  LASA  and  ALQ  which  are  very  close  to  the 
NTS.  For  the  travel  paths  to  LASA  and  ALQ,  the  deviation  of 
the  group  velocities  from  the  rest  of  the  travel  paths  might  be 
partially  due  to  the  fact  that  the  epicentral  distance  is  short. 
For  such  a short  epicentral  distance,  the  surface  waves  might 
not  be  well  dispersed  and  the  calculated  group  velocities  are 
very  sensitive  to  the  accuracy  of  the  event  locations  and  the 
origin  times  and  to  the  analyst's  arrival  time  picks.  In  addi- 
tion, for  these  travel  paths,  the  group  velocities  exhibit  with 
varying  degree  reverse  dispersion  (i,  e.  , group  velocity  de- 
creasing with  increasing  period)  at  periods  above  50  seconds. 

Referring  to  Figure  III-4,  which  gives  the  Love  wave  group 
velocities  for  the  travel  paths  from  NTS  to  the  observation  stations,  the 
above  remarks  for  the  Rayleigh  wave  group  velocities  along  the  same  travel 
paths  are  thought,  in  general,  to  be  applicable  here.  However,  it  is  noticed 
that  the  group  velocities  for  the  travel  paths  to  CTA,  KIP,  MAT,  and  GUM, 
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which  are  entirely  oceanic,  are  more  irregular  and  deviate  more  from  those 
for  the  typical  oceanic  path  as  compared  to  those  for  the  Rayleigh  waves  along 
the  same  travel  paths.  Considering  that  the  observed  Rayleigh  wave  group  ve- 
locities can  be  well  explained  with  a typical  ocean  path  Rayleigh  wave  group 
velocity  curve,  the  large  deviations  of  observed  Love  wave  group  velocities 
from  those  for  the  typical  oceanic  path  seem  surprising.  Although  the  phys- 
ical irregularities  involved  in  these  travel  paths  are  not  known,  these  ob- 
served irregularities  in  the  obtained  Love  wave  group  velocities  might  be  due 
to  the  inherent  nature  of  the  Love  wave  propagation  which  is  more  sensitive 
to  travel  path  irregularities  and  is  more  susceptible  to  multipath  effects  than 
is  the  Rayleigh  wave. 

Figures  III- 5 and  III-6  show  measured  Rayleigh  and  Love  wave 
group  velocities  from  EKA  to  the  observation  stations.  As  compared  to  trav- 
el paths  from  NIS  to  the  observation  stations,  there  are  fewer  events  used  to 
form  the  average  group  velocities  for  these  travel  paths.  In  fact,  the  Rayleigh 
wave  group  velocities  for  half  of  the  event- station  paths  shown  in  Figure  UI-5 
are  obtained  from  just  one  EKZ  event,  and  the  Love  wave  group  velocities  as 
shown  in  Figure  III- 6 are  obtained  only  for  five  event- station  paths.  That  is 
due  to  the  fact  that  there  are  fewer  events  and  observation  stations  available 
for  EKZ.  Also,  Love  wave  recordings  for  the  EKZ  events,  in  general,  are 
scarce.  Nevertheless,  the  Rayleigh  wave  group  velocities  shown  in  Figure 
111-5  are  thought  to  be  well  behaved,  even  for  those  event- station  paths  where 
the  group  velocities  are  obtained  from  just  one  EKZ  event.  It  is  obvious,  how- 
ever, that  the  travel  path  from  EKZ  to  CHG  shows  very  low  group  velocities. 
These  low  group  velocities  are  consistent  with  the  previous  observations  at 
CHG  for  the  central  Eurasian  earthquakes  and  the  eastern  Kazakh  explosions 
(Turnbull, et  al,  , 1973).  The  travel  paths  from  central  Eurasia  to  CHG  en- 
counter fairly  complicated  geological  features  in  central  Eurasia:  the  ab- 
normally thick  crust  (about  60  km  thickness)  of  the  Tibet  platform,  and  the 
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Himalayan  fold  system.  Thus,  these  low  group  velocities  are  thought  to  be 
normal  for  the  central  Eurasia  to  CHG  travel  paths. 

The  average  group  velocities,  as  shown  in  Figures  III-7  and 
ni-8,  for  the  travel  paths  from  the  north  Caspian  Sea  region  to  the  observa- 
tion stations  are  determined  from  just  three  PNE  events.  Except  for  the  trav- 
el paths  to  NORSAR  and  LASA,  the  group  velocities  for  the  other  four  travel 
paths  as  indicated  in  Figure  UI-7  represent  the  results  from  just  one  PNE 
event.  However,  judging  from  the  fairly  smooth  variation  of  these  group 
velocities,  similar  results  probably  can  be  expected  for  the  same  travel  paths 
from  examining  more  events  in  this  region. 

B.  OBSERVED  SURFACE  WAVE  AMPLITUDE  SPECTRA 
1.  Spectral  Contents  and  Spectral  Ratios 

The  observed  surface  wave  amplitude  spectra  at  all  available 
observation  stations  for  all  selected  events  listed  in  Table  II- 1 have  been  ob- 
tained in  the  manner  discussed  in  Subsection  II-C  and  are  displayed  in  Figures 
III- 9 through  III- 35.  For  each  event,  the  observed  spectra  have  been  equalized 
at  all  observation  stations  in  order  to  make  direct  comparison  possible.  The 
spectral  equalization  means  that  raw  observed  surface  wave  spectra  have 
been  corrected  for  the  station's  instrument  response,  geometric  spreading, 
and  travel  path  attenuation.  Hence,  the  spectra  shown  in  Figures  III- 9 through 
III- 35  are  the  spectra  as  they  would  be  observed  at  the  source.  For  the  trav- 
el path  correction,  Tryggvason's  attenuation  correction  (1965)  has  been  em- 
ployed for  all  event- station  paths,  except  for  the  travel  paths  from  EKZ  and 
PNE  to  CHC  where  the  attenuation  coefficients  calculated  specifically  for 
these  travel  paths  (Turnbull, et  al.  , 1973)  have  been  used.  The  reason  for  us- 
ing Tryggvason's  correction  will  be  given  later  when  we  discuss  the  travel 
path  attenuation  in  Section  IV.  In  each  figure,  available  observed  Rayleigh 
and  Love  wave  spectra  are  displayed  in  part  (a),  side-by-side  for  each 
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FIGURE  III- 9 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/926/74 
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FIGURE  III-9 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/926/74 
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FIGURE  III- 10 

OBSERVED  SURFACE  WAVE  SPECTRA;  NTS/430/75 
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(b)  LQ/LR  Amplitude  Spectral  Ratios 
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FIGURE  III- 10 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/430/75 
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OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/514/75 
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(a)  Amplitude  Spectra 
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FIGURE  III- 11 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/514/75 

(PAGE  2 OF  5) 
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FIGURE  IIMl 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/514/75 
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FIGURE  m-1 1 

OBSERVED  SURFACE  WAVE  SPECTRA;  NTS/514/75 
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FIGURE  III- 12 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTl/603/75 
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FIGURE  III-12 

OBSERVED  SURFACE  WAVE  SPECTRA: 
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FIGURE  III- 12 

OBSERVED  SURFACF  WAVE  SPECTRA:  NTl/603/75 
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FIGURE  III- 13 

OUSERVED  SURFACE  WAVE  SPECTRA:  NT2/b03/75 
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FIGURE  III- 14 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/619/75 
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(a)  Amplitude  Spectra 


FIGURE  III- 14 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/619/75 
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(a)  Amplitude  Spectra 
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FIGURE  III- 14 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/619/75 
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(b)  LQ/LR  Amplitude  Spectral  Ratio 
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FIGURE  m-14 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/619/75 
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(b)  LQ/LR  Amplitude  Spectral  Ratio 
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FIGURE  III- 14 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/619/75 
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(a)  Amplitude  Spectra 
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FIGURE  III- 15 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/626/75 
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(b)  LQ/LR  Amplitude  Spectral  Ratio 


■o 


IQ.Q  25.Q  40.0 

T C5EC] 

FIGURE  III- 15 
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FIGURE  III- 15 

OBSERVED  SURFACE  WAVE  SPECTRA;  NTS/626/75 

(PAGE  4 OF  4) 


(a)  Amplitude  Spectra 


10<0  25.  Q 40.  Q 10.0  ^S.O  40.0 


T C5ECJ  T CSEC] 

FIGURE  III- 16 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/1028/5 
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(a)  Amplitude  Spectra 
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FIGURE  UI-16 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/1028/5 
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FIGURE  III- 16 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/1028/5 
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(b)  LQ/LR  Amplitude  Spectral  Ratio 


FIGURE  III- 16 

OBSERVED  SURFACE  WAVE  SPECTRA;  NTS/1028/5 
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(b)  LQ/LR  Amplitude  Spectral  Ratio 


T C5ECD 
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FIGURE  ni-17 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/103/76 
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(b)  LQ/LR  Amplitude  Spectral  Ratio 
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FIGURE  III- 17 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/103/76 
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FIGURE  III- 18 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTl/204/76 
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(b)  LQ/LR  Amplitude  Spectral  Ratio 
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FIGURE  m-18 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTI/204/76 
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FIGURE  m-19 

OBSERVED  SURFACE  WAVE  SPECTRA:  NT2/204/76 
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FIGURE  III- 20 

I OBSERVED  SURFACE  WAVE  SPECTRA;  NTS/212/76 
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FIGURE  III- 20 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/212/76 

(PAGE  2 OF  4) 
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(b)  LQ/LR  Amplitude  Spectral  Ratio 
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FIGURE  III- 20 

OBSERVED  SURFACE  WAVE  SPECTRA;  NTS/212/76 
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(a)  Amplitude  Spectra 
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(a)  Amplitude  Spectra 
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FIGURE  III- 21 

OBSERVED  SURFACE  WAVE  SPECTRA;  NTS/309/76 

(PAGE  3 OF  5) 
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FIGURE  III- 21 


OBSERVED  SURFACE  WAVE  SPECTRA;  NTS/309/76 
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FIGURE  m-21 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/309/76 
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FIGURE  III- 22 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/314/76 
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FIGURE  III- 22 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/314/76 
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FIGURE  m-22 

OBSERVED  SURFACE  WAVE  SPECTRA;  NTS/314/76 

(PAGE  5 OF  6) 
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(b)  LQ/LR  Amplitude  Spectral  Ratio 
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FIGURE  III- 22 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTS/314/76 
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FIGURE  III- 23 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTl/317/76 
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FIGURE  III- 23 

OBSERVED  SURFACE  WAVE  SPECTRA:  NTl/317/76 
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FIGURE  m-24 

OBSERVED  SURFACE  WAVE  SPECTRA:  NT2/317/76 
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FIGURE  III- 24 

OBSERVED  SURFACE  WAVE  SPECTRA:  NT2/317/76 
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FIGURE  III- 24 

OBSERVED  SURFACE  WAVE  SPECTRA;  NT2/317/76 
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(b)  LQ/LR  Amplitude  Spectral  Ratio 
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FIGURE  III- 24 

OBSERVED  SURFACE  WAVE  SPECTRA:  NT2/317/76 
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(b)  LQ/LR  Amplitude  Spectral  Ratio 
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FIGURE  III- 24 

OBSERVED  SURFACE  WAVE  SPECTRA:  NT2/317/76 
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FIGURE  111-26 

OBSERVED  SURFACE  WAVE  SPECTRA;  EKZ/1026/3 
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FIGURE  III- 26 

OBSERVED  SURFACE  WAVE  SPECTRA: 
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FIGURE  III- 27 

OBSERVED  SURFACE  WAVE  SPECTRA:  EKZ/1214/3 
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(b)  LQ/LR  Amplitude  Spectral  Ratio 
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FIGURE  III- 27 

OBSERVED  SURFACE  WAVE  SPECTRA:  EKZ/1214/3 
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FIGURE  III-28 

OBSERVED  SURFACE  WAVE  SPECTRA:  EKZ/220/75 
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(b)  LQ/JLR  Amplitude  Spectral  Ratio 
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FIGURE  III- 28 

OBSERVED  SURFACE  WAVE  SPECTRA;  EKZ/220/75 
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(b)  LQ/LR  Amplitude  Spectral  Ratio 
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FIGURE  III- 29 

OBSERVED  SURFACE  WAVE  SPECTRA:  EKZ/311/75 
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FIGURE  UI-30 

OBSERVED  SURFACE  WAVE  SPECTRA;  EKZ/427/75 
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FIGURE  III- 32 

OBSERVED  SURFACE  WAVE  SPECTRA:  EKZ/704/76 
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FIGURE  III- 32 

OBSERVED  SURFACE  WAVE  SPECTRA;  EKZ/704/76 
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OBSERVED  SURFACE  WAVE  SPECTRA:  EKZ/704/76 
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j OBSERVED  SURFACE  WAVE  SPECTRA;  PNE/J222/I 

(PAGE  1 OF  2) 

I 


40.0 


I 


111-113 


(b)  LQ/LR  Amplitude  Spectral  Ratio 
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FIGURE  III- 33 

OBSERVED  SURFACE  WAVE  SPECTRA;  PNE/1222/1 
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FIGURE  III- 34 

OBSERVED  SURFACE  WAVE  SPECTRA:  PNE/820/72 
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observation  station  of  a given  event.  The  Love  wave  to  Rayleigh  wave  ampli- 
tude ratios  (LQ/LR),  if  available,  are  given  in  part  (b).  It  is  noticed  that 
these  observed  spectra  are  represented  by  spectral  values  at  nine  discrete 
periods,  namely  1 0 to  50  seconds  at  5 second  increments,  and  can  be  con- 
sidered as  the  straight-line-interpolated  version  of  the  continuous  spectra 
for  the  period  range  from  10  to  50  seconds. 

a.  The  NTS  events 

Figures  III-9  through  III-24  show  the  observed  surface  wave 
amplitude  spectra  and  LQ/LR  spectral  ratios  for  the  selected  NTS  events. 

The  results  of  examining  these  spectra  and  LQ/LR  spectral  ratios  can  be 
summarized  as  follows: 

• For  all  events,  most  of  the  observation  stations  indicate  Love 
wave  motion,  although,  theoretically,  an  explosive  source  in 
a horizontally  layered  homogeneous  isotropic  medium  should 
not  generate  Love  waves.  These  observations  of  the  Love 
waves  from  an  explosion  event  can  be  best  explained  by  the 
tectonic  strain  release  associated  with  the  explosion  (Toksdz, 
et  al,  , 1965).  However,  we  have  not  observed  the  situation 
that  the  Love  wave  amplitude  is  always  smaller  than  the  cor- 
responding Rayleigh  wave  amplitude  at  a given  observation  sta- 
tion as  Lambert  et  al.  (1972)  observed.  This  is  thought  to  be 
reasonah;le  since,  at  some  observation  stations,  Rayleigh  waves 
generated  from  a point  explosion  source  might  destructively 
interfere  with  the  Rayleigh  waves  radiated  from  the  tectonic 
strain  release,  while  the  Love  wave  radiated  from  the  tectonic 
strain  release  might  be  at  the  maximum  of  its  radiation  pat- 
tern. 

• The  shapes  of  the  observed  spectra  vary  from  event  to  event 
for  any  given  observation  station.  It  is  not  apparent  from 
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visual  inspection  whether  these  events  have  similar  source 
mechanisms.  However,  there  is  a common  feature  in  these 
spectra  with  a few  exceptions.  The  common  feature  is  that  the 
average  spectral  amplitude  is  higher  at  short  periods  and  de- 
creases toward  longer  periods.  From  knowledge  of  the  theo- 
retical spectral  variation  with  respect  to  focal  depth  for  the 
explosive  and  double-couple  source,  this  kind  of  spectral  be- 
havior is  a definite  indication  of  very  shallow  focal  depth.  This 
conjecture  of  very  shallow  focal  depth,  that  will  be  shown  later 
in  the  source  parameter  estimations  of  these  events  in  Section 
V,  justifies  the  assumption  made  by  Lambert  et  al.  (1972), 
Toksdz  et  al.  (1964,  1965)  and  Mitchell  (1975)  that  the  double- 
couple source  modeling  the  explosion-associated  tectonic  strain 
release  is  at  the  shot  depth  of  the  explosion. 

For  any  given  event,  the  level  of  the  observed  spectra  varies 
from  observation  station  to  observation  station.  This  azimuthal 
variation  of  the  observed  spectra  can  be  seen  more  clearly 
when  the  observed  radiation  patterns  are  plotted  in  Subsection 
V-C  for  several  events  with  good  azimuthal  coverage.  This 
variation  of  the  spectral  level  among  the  various  observation 
stations  indicates  that  surface  wave  radiation  from  these  events 
is  not  isotropic,  i.  e.  , the  radiation  pattern  is  non-circular.  In 
addition  to  the  observation  of  Love  waves,  this  is  another  in- 
dication of  possible  tectonic  strain  release  associated  with 
these  events.  It  is  also  noticed  that,  for  any  given  observa- 
tion station,  the  spectral  level  of  the  observed  spectra  varies 
from  event  to  event.  This  is  definitely,  although  not  entirely, 
due  to  the  fact  that  these  events  have  different  sizes,  as  in- 
dicated by  their  m^.  However,  this  also  might  be  an  indica- 
tion that  double-couple  sources  modeling  the  possible  tectonic 
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strain  release  associated  with  these  explosion  events  might 

not  all  have  the  same  strike  directions.  This  conjecture  of 

different  strike  directions,  that  will  be  shown  later  in  the 

source  parameter  estimations  of  these  events  in  Section  V, 

does  not  quite  agree  with  the  assumption  made  by  some  au- 

I 

thors  (Lambert,  et  al.  , 1972;  Toksoz,  et  al,  , 1964,  1965) 
that  the  double-couple  source  modeling  of  the  explosion- 
associated  tectonic  strain  release  has  the  same  strike  direc- 
tion as  that  of  the  local  earthquakes  in  the  same  region, 

• For  some  events,  it  is  noticed  that  the  spectral  level  of  the 
observed  spectra  at  the  SDCS  stations  is  two  to  three  times 
higher  than  that  at  other  stations.  There  is  no  known  physical 
explanation  for  this  difference,  although  it  is  suspected  that 
this  relatively  high  spectral  level  is  spurious  and  probably  due 
to  the  miscalibration  of  the  instruments  in  the  field  recordings 
or  in  the  analog-to-digital  (A/D)  conversion.  The  reason  for 
this  suspicion  is  based  on  the  examination  of  the  signal  cali- 
bration value  (m^/cc)  which  shows  abnormally  large  variation 
from  event  to  event  and  from  trace  to  trace,  ranging  from  1 
mfi/cc  to  20  mu/cc.  However,  the  spectral  shapes  of  the  ob- 
served spectra  at  the  SDCS  stations  seem  comparable  to  those 
at  other  stations.  Therefore,  in  the  source  parameter  esti- 
mations of  these  events  in  Section  V,  the  observed  spectra 
from  the  SDCS  stations  are  only  used  for  those  events  for 
which  no  other  data  are  available. 

• The  LQ/LR  spectral  ratios,  as  shown  in  part  (b)  of  Figures 
III-9  through  III-24,  range  from  0,  if  we  assign  LQ/LR  = 0 
for  those  observation  stations  which  do  not  have  useable  Love 
wave  recordings,  to  about  10,  taking  into  account  all  events. 
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stations,  and  periods.  These  LQ/LR  spectral  ratios  are  com- 
patible with  the  average  fundamental  mode  LQ/LR  spectral 
ratios  for  the  shallow  earthquakes  theoretically  modeled  by 
the  double-couple  sources  with  all  possible  source  mechanisms, 
i.  e.  , dip  angle  ranging  from  1 0 to  90  degrees,  slip  angle  rang- 
ing from  -90  to  +90  degrees,  and  strike  angle  ranging  from  0 
to  180  degrees  (Turnbull,  1976).  This  does  not  quite  agree 
with  the  observations  made  by  Lambert  et  al.  (1972)  that 
LQ/LR  spectral  ratios  for  explosions  are  less  than  one  and 
are  much  less  than  those  for  earthquakes.  Also,  it  is  noticed 
that,  on  the  average,  the  LQ/LR  spectral  ratios  increase,  al- 
though not  significantly,  from  10  seconds  toward  50  seconds. 
This  is  somewhat  different  from  the  theoretical  results  for 
earthquakes  given  by  Turnbull  (1976)  where  the  average  funda- 
mental mode  LQ/LR  spectral  ratios  increase  very  noticeably 
from  10  seconds  toward  50  seconds.  However,  this  difference 
is  not  thought  significant  enough  to  be  used  as  the  discrimina- 
tion feature  between  the  earthquake  and  the  explosion  as  im- 
plicitly suggested  by  Lambert  et  al.  (1972). 

• In  order  to  make  possible  comparison  of  the  LQ/LR  spectral 

ratios  among  the  selected  events,  it  is  helpful  to  tabulate  the 
values  of  the  LQ/LR  spectral  ratios  for  all  selected  events 
observed  at  the  same  station.  Referring  to  Table  II- 1,  it  is 
noticed  that  stations  KON  and  ALFA  have  both  Rayleigh  and 
Love  wave  recordings  for  most  of  the  selected  events.  There- 
fore, ^n  Table  III-l,  the  LQ/LR  spectral  ratios  observed  at 
KON  and  ALFA  are  given  with  the  events  ordered  in  decreas- 
ing event  m^.  These  values  in  Table  III-l  are  thought  to  be 
representative  for  the  LQ/LR  spectral  ratios  of  the  selected 
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events.  Theoretically,  if  two  events  have  the  same  or  very 
similar  focal  depth,  dip,  slip  and  strike  angle,  and  F value, 
their  LQ/LR  spectral  ratios  observed  at  the  same  observation 
station  should  be  very  similar.  Referring  to  Table  III- 1 , al- 
though no  exactly  identical  LQ/LR  spectral  ratios  are  observed 
among  these  events,  it  is  noticed  that  the  events  NTS/514/75, 
NTl/603/75,  and  NTS/619/75  have  similar  LQ/LR  spectral 
ratios.  Also,  the  events  NTS/1028/5  and  NTS/626/75  have 
similar  LQ/LR  spectral  ratios.  There  is  no  apparent  cor- 
relation between  the  event  m^  and  the  LQ/LR  spectral  ratio. 

This  is  expected,  since  the  difference  in  the  event  m should 

b 

disappear  in  the  process  of  taking  the  Love  wave  to  the  Ray- 
leigh wave  amplitude  spectral  ratio, 

b.  The  EKZ  events 

Figures  III-25  through  III-32  show  the  observed  surface  wave 
amplitude  spectra  and  the  LQ/LR  spectral  ratios  for  the  selected  EKZ  events. 
The  spectral  quality  of  these  EKZ  events  is  not  as  good  as  that  of  the  NTS 
events.  There  are  fewer  stations  where  good  surface  wave  recordings  are 
available.  An  exception  to  this  is  the  event  EKZ/704/76,  Also,  Love  wave 
recordings  are  not  as  frequently  available  as  for  NTS  events.  Nevertheless, 
for  all  selected  EKZ  events.  Love  wave  recordings  do  exist  at  least  for  one 
observation  station.  Whenever  they  are  available,  their  amplitude  spectia 
seem  compatible  with  the  corresponding  Rayleigh  wave  amplitude  spectra. 

In  general,  results  similar  to  those  stated  previously  for  the 
NTS  events  about  the  variation  of  the  spectral  shape  and  the  spectral  level 
are  found  for  the  selected  EKZ  events.  Again,  as  for  the  selected  NTS  events, 
the  shape  of  the  observed  amplitude  spectra  varies  from  event  to  event  for 
any  given  observation  station,  yet  the  spectral  shape  indicates  a very  shallow 
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event  depth.  Although  Love  wave  motions  are  observed  only  at  a few  observa- 
tion stations,  variation  of  the  spectral  level  of  observed  Rayleigh  wave  ampli- 
tude spectra  does  exist  among  the  observation  stations  for  a given  event.  This 
azimuthal  variation  of  the  observed  Rayleigh  wave  amplitude  suggests  that  the 
selected  EKZ  events  do  contain  some  double-couple  component. 

The  LQ/LR  spectral  ratios,  as  shown  in  part  (b)  of  Figures 
111-25  through  III- 32,  range  from  0 to  5.  In  general,  these  LQ/LR  spectral 
ratios  seem  slightly  lower  than  those  for  the  NTS  events;  however,  the  dif- 
ference is  not  thought  to  be  significant  enough  to  make  a distinction  between 
the  explosions  at  EKZ  and  NTS.  The  typical  values  of  the  LQ/LR  spectral 
ratios  observed  at  stations  ALFA  and  NORS  are  given  in  Table  III- 2.  As 
with  the  NTS  events,  no  obvious  correlations  are  observed  between  these 
values  and  the  event  m^.  No  noticeable  similarity  in  the  LQ/LR  spectral 
ratios  among  the  events  is  observed. 

c.  The  PNE  events 

Figures  III-33  through  III-35  show  the  observed  surface  wave 
amplitude  spectra  and  the  LQ/LR  spectral  ratios  for  the  selected  PNE  events. 
Also  the  typical  values  of  the  LQ/LR  spectral  ratios  observed  at  stations 
NORS  and  1.ASA  are  shown  in  Table  III- 3.  Love  wave  motion  is  observed  at 
most  of  the  stations  for  these  three  PNE  events.  As  observed  in  the  NTS 
and  the  EKZ  events,  the  general  shapes  of  the  observed  amplitude  spectra 
suggest  that  these  events  took  place  at  very  shallow  depths.  Although  there 
are  only  two  observation  stations  involved  in  the  events  PNE/1222/1  and 
PNE/820/72,  variation  in  the  observed  surface  wave  amplitude  spectra  is 
noticeable  among  observation  stations  for  all  three  events.  This  anisotrop- 
ic radiation  of  the  surface  waves  together  with  the  observation  of  the  Love 
wave  motion  is  again  indicative  of  tectonic  strain  release  associated  with 
these  explosions.  The  range  of  the  variation  of  the  LQ/LR  spectral  ratios 
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TYPICAL  LQ/LR  SPECTRAL  RATIOS  FOR  THE  SELECTED  PNE  EVENTS 


Event  1.  D. 

Period  (sec) 

50 

45 

40 

35 

30 

25 

20 

15 

10 

Station:  NORS 

i 

PNE/729/76 

5.  9 

0,  92 

0.  40 

0.  27 

0.  24 

0.  29 

0.  19 

0.  18 

0.  40 

PNE/820/72 

5.7 

0.  73 

0.  55 

0.  36 

0.  33 

0.  29 

0,  32 

0.  35 

0.  26 

0.  71 

Station:  LASA 

PNE/1222/1 

6.  0 

3.  02 

2.  10 

4.  83 

0.  84 

0.  39 

0.  40 

0.  57 

0.  43 

1.  20 

PNE/729/76 

5.  9 

2.  55 

2.  84 

2.  67 

0.  79 

0.  37 

0,  28 

0.46 

0.  83 

2.  05 

PNE/820/72 

5,  7 

2.  51 

2.  69 

2.  69 

1. 05 

0,  31 



0.  41 

3.  09 
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is  very  much  the  same  as  that  for  the  EKZ  events,  from  0 to  about  5.  Refer- 
ring to  Table  ItI-3,  the  LQ/LR  spectral  ratios  for  the  event  PNE/820/72  are 
quite  similar  to  those  for  the  event  PNE/729/76,  indicating  that  these  two 
events  very  likely  have  similar  source  parameters  in  terms  of  depth,  dip, 
slip  and  strike  angle,  and  the  F value. 

2.  Relative  Love  Wave  Excitation  from  a Suite  of  NTS  Events 


As  discussed  in  Subsection  lII-B-l-a,  Love  wave  motion  is 

observed  for  all  selected  NTS  events  at  most  observation  stations.  However, 

at  any  given  observation  station,  no  apparent  correlation  between  the  level  of 

the  observed  Love  wave  amplitude  spectra  and  the  event  m is  observed. 

b 

This  is  thought  to  be  reasonable,  since  at  any  given  observation  station  the 
observed  spectral  level  of  the  Love  wave  amplitude  spectra  not  only  depends 
on  the  event  m^  but  also  depends  (even  more  heavily)  on  the  F value,  which 
is  the  indication  of  the  amount  of  the  tectonic  strain  release  associated  with 
the  explosion,  and  the  strike  direction  of  the  double-couple  source  which 
models  the  explosion-associated  tectonic  strain  release.  For  example,  even 
assuming  that  two  explosion  events  have  the  same  amount  of  tectonic  strain 
release,  the  Love  wave  amplitude  spectral  levels  for  these  two  events  ob- 
served at  the  same  station  can  be  very  different  just  due  to  the  different  ori- 
entations of  the  Love  w’ave  radiation  patterns  caused  by  the  different  strike 
directions. 


To  ex  imine  the  relative  Love  wave  excitation  from  a suite  of 
NTS  events,  two  approaches  are  taken.  In  the  first  approach,  one  NTS  event 
is  taken  as  a reference  event  and  the  relative  Love  wave  excitation  is  deter- 
mined from  the  ratios  of  the  observed  Love  wave  amplitude  spectra  of  a given 
event  to  those  of  the  reference  event.  These  ratios  are  calculated  for  all 
available  common  observation  stations  where  they  are  compared  among  select- 
ed NTS  events.  In  the  second  approach,  the  seismic  moment  of  the  tectonic 
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strain  release  associated  with  the  explosion  is  estimated  for  each  NTS  event 
and  these  estimated  seismic  moments  are  then  compared  among  the  selected 
events. 

For  the  first  approach,  the  event  NTS/619/75,  which  has  ob- 
served Love  wave  amplitude  spectra  at  all  six  observation  stations,  is  chosen 
as  the  reference  event.  The  Love  wave  amplitude  ratios  of  the  other  NTS 
events  to  the  reference  event  NTS/619/75  are  shown  in  Table  III-4  for  the  ob- 
servation stations  ALFA  and  KON.  In  Table  III-4,  the  events  are  arranged 

in  the  order  of  decreasing  m,  . Theoretically,  if  these  events  have  identical 

b 

or  very  similar  source  parameters,  then  their  Love  wave  ratios  should  be 
very  similar  at  any  observation  station  and  these  ratios  can  be  conveniently 
used  as  a measure  of  the  relative  Love  wave  excitations  (Lambert,  et  al.  , 
1972).  However,  without  knowledge  of  the  source  parameters  of  these  events, 
the  direct  comparisons  of  these  Love  wave  ratios  can  be  very  misleading. 
Referring  to  Table  III-4,  it  is  obvious  that  quite  different  conclusions  about 
the  relative  Love  wave  excitation  from  these  events  can  be  drawn  by  examin- 
ing the  Love  wave  ratios  observed  at  station  KON  or  by  examining  those  ob- 
served at  station  ALFA,  For  example,  for  event  NTS/430/75,  the  Love  wave 
ratios  observed  at  ALFA  indicate  that  this  event  should  have  much  less  Love 
wave  excitation  than  event  NTS/619/75;  while  the  Love  wave  ratios  observed 
at  station  KON  suggest  that  the  Love  wave  excitation  from  event  NTS/430/75 
should  be  about  1.  5 times  larger  than  that  from  event  NTS/619/75,  These 
differences  in  Love  wave  ratios  for  the  same  event  observed  at  different  ob- 
servation stations  can  not  be  explained  by  the  m^  differences  between  event 
NTS/430/75  and  the  reference  event  NTS/619/75.  Moreover,  the  Love  wave 
ratios  observed  at  the  same  observation  stations  for  events  with  comparable 

m , such  as  the  events  NTS/626/75  and  NTS/103/76,  show  drastic  differences, 
b 

No  apparent  correlation  between  the  Love  wave  ratios  and  the  event  m^  are 
observed,  i.  e.  , the  larger  m,  does  not  necessarily  imply  the  larger  Love 
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LOVE  WAVE  AMPLITUDE  RATIOS  OF  THE  SELECTED  NTS  EVENTS: 
REFERENCE  EVENT  NTS/619/75  m.  = 6.  1 

D 


Event  1.  D. 

m 

b 

Period  (sec) 

mm 

40 

35 

30 

25 

20 

15 

10 

Station:  KCN 

NTS/212/76 

6,  3 

2.  59 

3.  40 

2.  85 

2.  96 

3.  84 

3.  62 

2.  58 

2.  25 

1.  20 

NTS/626/75 

6.  2 

1. 76 

1. 87 

1. 70 

1. 69 

1. 78 

1. 67 

1.  55 

1.  55 

0.  57 

NTl  /31  7/76 

6.  1 

0.  75 

0.  84 

0.  65 

0.  60 

0.  50 

0.  16 

0.  23 

0.  47 

0.  09 

NTS/514/75 

6.  0 

0.  56 

0,  62 

0.  63 

0.  68 

1. 05 

0.  96 

0.  96 

0.  52 

0.  36 

NTl /603/75 

5.  9 

0.  36 

0.  44 

0,  39 

0.48 

0.  66 

0.  58 

0.  52 

0.  31 

0.  09 

NT2/31  7/76 

5.  8 

0.  32 

0.  37 

0.  35 

0.  24 

0.  22 

0.  29 

0.  41 

0.  47 

0.  16 

NT2/603/76 

5.  7 

0.  20 

0.  22 

0.  21 

0.  21 

0.  27 

0.  29 

0.  42 

0.  31 

0.  14 

NTS/430/75 

5.  2 

1. 09 

1.  33 

1.  21 

1.  56 

2.  03 

1.  50 

2.  16 

1. 46 

1.  15 

Station:  ALPA 

NTS/1028/5 

6.  4 

5.  02 

5.  38 

5.  02 

4.  57 

4.  37 

4.  23 

4.  51 

2.  27 

2.  71 

NTS/626/75 

6.  2 

2.  87 

2,  88 

2,  72 

2.  50 

2.  43 

2.  41 

2.  58 

2.  08 

1.13 

NTS/103/76 

6.  2 

9.  64 

9.  51 

8.  86 

8.  33 

7.  88 

7.  29 

8.  13 

4.  92 

5.  11 

NTS/514/75 

0.  69 

0.  79 

0.  76 

0,  74 

0.  71 

0.  64 

0.  69 

0.  28 

0.  62 

NTl  /603/75 

0.  23 

0.  27 

0.  28 

0.  26 

0.  29 

0.  27 

0.  21 

0.  29 

NT2/603/75 

0.  12 

0.  09 

0.  10 

0.  09 

0.  08 

0.  13 

0.  09 

0.  12 

NTS/926/74 

0.  05 

0.  03 

0.  03 

0.  02 

0.  03 

0.  06 

0.  04 

0.  12 

NTS/430/75 

5.  2 

0.  18 

0.  07 

0.  05 

0.  03 

0.  03 

0.  03 

0.  04 

0.  04 

wave  ratios.  It  is  believed  that  these  large  differences  in  the  observed  Love 
wave  ratios  are  probably  caused  by  the  differences  in  the  F values  and  the 
source  parameters  among  the  selected  NTS  events.  Therefore,  it  seems  not 
too  plausible  to  obtain  reliable  information  about  the  relative  Love  wave  ex- 
citation by  direct  comparison  of  amplitude  spectra  or  amplitude  ratios  ob- 
served at  a number  of  common  observation  stations. 

For  the  second  approach,  Table  III- 5 presents  the  estimated 

seismic  moments  of  the  tectonic  strain  releases  associated  with  the  selected 

NTS  events,  their  locations,  and  event  m . The  details  of  how  these  seismic 

b 

moments  were  estimated  will  be  given  in  Section  V where  the  source  param- 
eter estimations  of  all  selected  NTS,  EKZ,  and  PNE  events  are  discussed. 

In  Table  III-5,  the  events  are  ordered  in  decreasing  seismic  moment.  Thus. 
Table  III- 5 actually  gives  the  rankings  of  the  selected  NTS  events  in  terms  of 
their  associated  tectonic  strain  release.  The  seismic  moment  is  estimated 
from  the  observed  surface  wave  amplitude  spectra  recorded  at  all  observa- 
tion stations.  Since  it  is  a constant  for  a given  event,  the  seismic  moment, 
which  measures  tectonic  strain  release  associated  with  an  explosion,  is 
thought  to  be  a more  stable  and  useful  quantity  than  direct  observations  of 
surface  wave  amplitude  spectra  or  Love  wave  spectral  ratios.  Referring  to 
Table  III- 5,  for  the  NTS  PAHUTE  events,  there  is  an  apparent  trend  between 
the  estimated  seismic  moments  of  the  explosion- as sociated  tectonic  strain 
releases  and  the  event  For  those  events  which  took  place  in  hard  rock, 

i.  e,  , shot  location  indicated  by  NTS  PAHUTE,  the  trend  indicates  that  the 
seismic  moment  of  the  explosion-associated  tectonic  strain  release  decreases 
with  decreasing  m^.  This  trend  is  thought  to  be  reasonable,  if  the  source  of 
the  tectonic  strain  energy  released  by  an  explosion  is  by  relaxation  of  the 
cavity  (Toksbz,  et  al.  , 1965;  Press  and  Archambeau,  1962).  One  event  which 

deviates  drastically  from  the  rest  of  the  events  is  NTS/430/75.  The  m of 

b 

this  event  is  5.  2 and  yet  the  estimated  seismic  moment  associated  with  this 
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TABLE  III-5 


ESTIMATED  SEISMIC  MOMENTS  OF  THE  TECTONIC  STRAIN  RELEASES 
ASSOCIATED  WITH  THE  SELECTED  NTS  EVENTS 


Event  1.  D. 

""b 

25 

M^  X 10  dyne-cm 

Location 

NTS/1028/5 

6.4 

0.  243 

NTS  PAHUTE 

NTS/626/75 

6.  2 

0.  127 

NTS  PAHUTE 

NTS/103/76 

6.  2 

0.  1 1 1 

NTS  PAHUTE 

NTS/212/7b 

6.  3 

0.  105 

NTS  PAHUTE 

NTS/314/76 

6.  3 

0. 0975 

NTS  PAHUTE 

NTl  /31  7/76 

6.  1 

0. 0428 

NTS  PAHUTE 

NTS/309/76 

6.  0 

0, 0332 

NTS  PAHUTE 

NTS/430/75 

5.  2 

0. 0212 

NTS 

NTS/619/75 

6.  1 

0.  0192 

NTS  PAHUTE 

NT2/317/76 

5.  8 

0.  0180 

NTS 

NTS/514/75 

6.  0 

0.  0170 

NTS  PAHUTE 

NTl  /603/75 

5.  9 

0. 00864 

NTS  PAHUTE 

NT2/603/75 

5.  7 

0. 00302 

NTS 

NTS/926/74 

5.  6 

0. 00163 

NTS 

*This  information  is  obtained  from  the  announced  United 
States  Nuclear  Test  Statistics  given  by  Nevada  Opera- 
tion Office,  Energy  Research  and  Administration 
(30  June  1 976) 
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explosion  is  what  is  expected  of  the  events  with  around  6,  0.  The  real 

physical  reason  for  this  large  deviation  is  not  known  at  the  present  time.  One 
possible  explanation  for  this  situation  is  that  a real  earthquake  was  triggered 
by  this  explosion.  In  that  case,  the  size  of  the  triggered  earthquake  is  not 
directly  dependent  on  the  size  of  the  explosion  (Brune  and  Pomeroy,  1963; 
Aki,  1964). 
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SECTION  IV 

TRAVEL  PATH  ATTENUATION 


. Several  factors  affect  the  spectral  level  and  shape  of  the  far- 
field  surface  wave  amplitude  spectra  of  a seismic  event.  These  factors  can 
be  conveniently  lumped  into  three  groups;  the  source  parameters  of  the 
event,  travel  path  effects,  and  the  instrument  response  of  the  observation 
station.  Therefore,  in  order  to  estimate  the  seismic  source  parameters  us- 
ing the  observed  far-field  surface  wave  spectra,  some  appropriate  corrections 
to  the  observed  spectra  must  be  made.  The  corrections  for  the  instrument 
response  and  geometric  spreading  are  trivial,  while  the  correction  for  the 
energy  attenuation,  which  is  the  other  part  of  the  travel  path  effect,  generally 
is  not  well  known  due  to  limited  knowledge  of  the  surface  wave  energy  dissipa- 
tion. 

The  travel  path  attenuation  correction  is  important  in  the  esti- 
mation of  the  seismic  source  parameters  by  spectral  fitting  of  observed  far- 
field  surface  wave  data.  The  level  of  the  attenuation  curve  can  affect  the 
estimation  of  the  seismic  moment,  and  the  general  shape  of  the  curve  can 
affect  the  estimation  of  the  focal  depth  and  the  source  mechanism.  So  far, 
when  no  reliable  knowledge  of  energy  attenuation  along  the  travel  paths  of 
interest  is  available,  Tryggva son's  energy  attenuation  curve  (Tryggvason, 

1965)  has  been  used.  The  general  shape  of  Tryggvason's  curve  is  believed 
to  be  reasonable,  although  the  level  of  the  curve  seems  a little  low  (Tryg- 
gvason, 1965;  Anderson,  1964).  Nevertheless,  in  the  past,  using  Tryggvason's 
attenuation  correction,  the  seismic  moment  estimates  of  numerous  Eurasian 
seismic  events  have  fitted  very  well  with  the  seismic  moment- versus -body wave 
magnitude  curve  determined  from  the  revised  w^-mod«l  given  by  Tsai  (1972b). 
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It  is  the  purpose  of  this  section  to  obtain  the  surface  wave  attenv'ation  curves 
for  several  travel  paths  of  interest  and  to  find  out  how  applicable  Tryf^gvason' s 
curve  is. 

Two  methods  to  determine  the  surface  wave  attenuation  coef- 
ficients will  be  derived  and  discussed  in  Subsection  A.  In  Subsection  B,  these 
methods  are  applied  to  the  observed  surface  wave  data  available  for  the  select- 
ed NTS  and  EKZ  events. 

A.  METHODS 

The  term  'energy  dissipation'  or  'energy  attenuation'  used 
throughout  this  report  is  chosen  to  indicate  all  energy  lost  along  the  travel 
path  including  absorption,  reflected  energy,  and  scattered  energy. 


I In  the  following  derivation,  it  is  assumed  that  there  is  no  energy 

\ transferred  from  one  frequency  to  another  along  the  travel  path  and  that  the 

I 

observed  amplitude  spectra  have  been  corrected  for  the  station  instrument 
response.  Two  methods  will  be  described  below:  the  two- station  method 
which  is  suitable  for  two  stations  along  the  same  event- station  travel  path, 
and  the  isotropic  - sou  rce  method  which  can  be  employed  for  the  explosion  data. 

1.  Two-Station  Method 

The  observed  surface  wave  amplitude  spectra  of  a given  seis- 
mic event  recorded  at  two  stations  can  be  related  as  follows: 


whe  re 


Aj(f) 


= R(0 


e ^ ^(IV-1) 


A.(f) 


is  the  observed  amplitude  spectrum  for  frequency 
f at  station  i. 
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1 
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k 


D. 

1 

4 

1 

O. 

1 

R(0j,02.f) 


is  the  amplitude  attenuation  coefficient  for  fre- 
quency f along  the  travel  path  i (great  circle 
path  between  the  seismic  event  and  station  i ), 
is  the  epicentral  distance  in  km  to  station  i , 
is  the  epicentral  distance  in  degrees  to  station  i, 
is  the  source  azimuthal  angle  to  station  i , and 
is  the  ratio  of  the  source  excitation  at  Oj  to 
that  at  O^,  for  frequency  f , 


In  equation  (IV- 1 ) the  term  RlO^.O^.f)  accounts  for  the  non- 
circular radiation  pattern  of  the  source.  For  an  isotropic  source,  R(0j,02,f) 
= 1.0.  For  an  anisotropic  source,  such  as  an  earthquake,  in  general 
R(Oj , O^,  f)  ¥ 1 . However,  in  the  special  case,  0^=0^.  R(0j,02,f)  = l re- 
gardless of  source  type.  That  is,  when  two  observation  stations  are  on  the 
same  great  circle  path  between  the  seismic  event  and  the  station  and  have  the 
same  azimuthal  angle,  equation  (IV- 1)  can  be  reduced  to 


Aj(f) 


(IV-2) 


where  k now  is  the  amplitude  attenuation  coefficient  for  frequency  f along 
A 

the  travel  path  between  stations  1 and  2.  Therefore,  equation  (IV-2)  can  be 
applied  to  the  surface  wave  data  recorded  at  two  stations  along  the  same 
travel  path  to  calculate  the  travel  path  attenuation  as  follows: 


k^(f) 


Aj(f) 


(IV- 3) 


Estimates  of  observational  errors  inherent  in  the  two- station 
method  might  be  inferred  in  the  following  way.  Hermann  (1974),  in  his  studies 
of  source  spectra  of  earthquakes  and  anelastic  attenuation  of  surface  waves. 
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found  that  the  source  spectra  could  be  estimated  within  a factor  of  1,  5 at  95% 
confidence  level  from  a good  data  set.  As  a consequence  of  this  error  esti- 
mate, one  can  state  that  in  such  a case  the  individual  spectral  amplitudes  at 
a recording  station  can  be  determined  within  a factor  of  1 . 5 accuracy,  at  the 
same  confidence  level.  Using  this  estimate  of  the  variation  in  the  observed 
spectra,  an  estimate  can  be  made  of  the  inherent  error  in  the  attenuation 
coefficient  determinations  by  the  two- station  method.  That  is,  in  the  worst 
possible  case  (the  amplitude  could  be  1.  5 times  larger  than  the  average 
and  A^  could  be  1.  5 times  smaller  than  the  average;  or  A^  could  be  1.  5 
times  smaller  and  A^  could  be  1 . 5 times  larger),  the  variation  of  the  cal- 
culated k . (f)  may  be  estimated  as: 

A 


2 f n(  1 . 
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2.  Isotropic-Source  Method 

When  data  from  more  than  two  stations  are  available,  two  ap- 
proaches other  than  the  two- station  m.ethod  can  be  used.  One  method  outlined 
by  Tsai  and  Aki  (1969)  assumes  that  the  source  characteristics  (focal  depth, 
source  mechanism,  and  average  structure  of  the  medium)  are  known;  while 
the  other  method  used  by  Tryggvason  (1965)  assumes  that  the  surface  wave 
radiation  is  isotropic.  There  are  advantages  and  disadvantages  in  using 
these  methods.  Tsai  and  Aki's  approach  can  take  into  account  possible  varia- 
tions in  the  observed  spectral  amplitudes  at  each  station  and,  hence,  can  be 
applied  to  the  observed  data  from  any  source  type  as  long  as  th<_  exact  source 
characteristics  are  known  a priori, which  is  the  shortcoming  of  this  method. 

On  the  other  hand,  Tryggvason's  approach  is  more  straight  forward  since  no 
prior  knowledge  about  the  source  characteristics  is  required.  However, 
Tryggvason's  method  can  only  be  applied  to  observed  data  from  an  isotropic 
source.  Also,  in  both  methods,  the  estimated  k^  value  represents  a region- 
al average  rather  than  a value  characteristic  of  a particular  travel  path  as  that 
obtained  through  the  two- station  method. 
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Here,  Tryggvason's  method  is  chosen  to  be  employed  by  the 
observed  surface  wave  data  from  the  selected  NTS  and  EKZ  events.  The 
main  reason  for  choosing  Tryggvason's  method  is  that  the  source  characteris- 
tics of  these  events  are  not  known,  and  explosion  events  are  closer  to  an  iso- 
tropic source  than  are  earthquake  events.  The  original  Tryggvason's  method 
only  applies  to  surface  wave  data  from  one  explosion  event.  In  the  following, 
Tryggvason's  method  will  be  extended  in  order  to  use  the  surface  wave  data 
from  many  explosion  events  from  the  same  area.  Also,  the  variance  of  the 

estimated  k ^ value  will  be  derived. 

A 

In  Tryggvason's  method,  the  surface  wave  amplitude  attenua- 
tion coefficient  can  be  estimated  from  the  following  equation  by  using  the  far- 
field  observed  surface  wave  amplitude  spectra  at  several  observation  stations 
for  one  explosion: 


fn  A.(f)  + |/n  sin4.  = C (f)-D.k  (f)  i=l,NSITE  (IV-5) 
1 1^1  /\ 

where  A.,  A.,  and  D.  were  previously  defined,  h (f)  is  now  the  average 
11  1 A 

amplitude  attenuation  coefficient  for  travel  paths  between  the  event  and  all 
NSITE  available  observation  stations  at  frequency  f . is  a constant  de- 

pending on  frequency  and  the  given  explosion  event.  In  equation  (IV-5),  k 

A 

and  are  two  unknown  parameters  to  be  estimated.  Hence,  when  NSITE 

is  greater  than  two,  k and  C can  be  estimated  by  linear  regression: 
i.  e.  , minimizing  the  following  quantity 


€ = 


N^TE  - 
1^1 


2.  (fnA.  + 4fnsin4.)  - (C-  - D 

1 ^ A 


(IV-6) 


The  estimated  will  be  as  follows: 

A 


k . (f) 
A 


NSITE  ^ + i sin/l.)  D.  ( /nA^  + ] f n sin 


? D.)^  - NSITE 
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where 


NS^E 

stands  for  A<  . It  can  be  easily  shown  that  equation  (IV- 7) 


1 i = l 

wnll  be  reduced  to  equation  (IV- 3)  when  NSITE  is  equal  to  two,  as  expected. 


By  regression  analysis,  the  standard  error  S of  the  estimated  values 

can  be  estimated  from  the  sample  as 


NSITE 
NSITE- 2 


^ D.  )^  - NSITE  ? 
11  11 
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We  have  discussed  the  method  for  obtaining  the  average  ampli- 
tude attenuation  coefficient  k together  with  its  error  estimate  S . This 

J\ 

process  can  be  performed  with  data  from  each  of  the  selected  events  in  the 

same  area  so  that  independent  estimates  of  k^  can  be  made  for  each  event. 

The  result  is  that  a mean  k value  can  be  obtained  together  with  its  error 

estimate.  Now  each  k . . has  been  determined  from  the  observed  data  at 

2 

NSITE.  stations  with  a corresponding  variance  S.  , where  j refers  to  the 
^ th 

particular  estimate  made  from  the  j explosion.  The  mean  k value  and 

2 ^ 

its  error  estimate  can  be  obtained  from  k.  . and  S.  as  follows  (Hoel,  1972; 

Aj  J 

Herrmann  and  Mitchell,  1975): 


NEXP 


j=l 


k , ./NEXP 
Aj 


(IV-9) 


NEXP 

S*^  = Y'  S'^/NEXP 

^ 


(IV-10) 


where  NEXI-’  is  the  number  of  available  explosions  in  tlie  same  area. 


B.  HE.St’I.TS 

I ()  calculate  the  surface  wave  amplitude  .ittenuation  coefficients 
for  the  travel  )ath  of  interest,  the  methods  discussed  in  tl  e previous  subsec- 
tion are  applied  to  the  Bayleigh  wave  data  from  two  groups  of  events,  namely 
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the  selected  NTS  and  EKZ  events  as  listed  in  Table  II- 1.  The  observed  Ray- 
leigh wave  amplitude  spectra  with  appropriate  instrument  response  correc- 
tion are  obtained  in  the  manner  described  in  Section  II.  In  the  following,  the 
results  will  be  presented  in  terms  of  two  groups  of  travel  paths,  i.  e.  , the 
travel  paths  from  the  NTS  to  the  observation  stations  available  to  the  select- 
ed NTS  events  and  the  travel  paths  from  EKZ  to  the  observation  stations  avail- 
able to  the  EKZ  events. 

1.  Travel  Paths  From  the  NTS  to  the  Observation  Stations 

Referring  to  Figure  III-l  which  shows  the  travel  paths  encounter- 
ed from  the  NTS  to  the  observation  stations,  these  travel  paths  can  be  logically 
divided  into  two  groups:  (1)  mostly  oceanic  p>aths  from  the  NTS  to  stations  in 
the  Pacific  Ocean  and  (2)  mostly  continental  or  mixed  (i.  e.  , partially  contin- 
ental and  partially  oceanic)  paths  from  the  NTS  to  the  rest  of  the  observation 
stations.  Therefore,  two  sets  of  Rayleigh  wave  attenuation  coefficients  are 
estimated  for  two  groups  oi  travel  paths.  The  travel  paths  in  the  first  group 
are  those  from  the  NTS  to  stations  CTA,  KIP,  and  MAT.  The  travel  paths 
in  the  second  group  include  those  from  the  NTS  to  stations  ALPA,  NORS, 

LASA,  CHG,  TLO,  KON,  OGD,  ALQ,  ZLP,  and  MAI.  The  travel  paths  from 
the  NTS  to  the  SDCS  Stations  are  purposely  excluded  due  to  the  non-uniform 
data  quality  at  those  stations  as  discussed  in  Section  111.  For  the  first  group 
of  travel  paths,  it  is  noticed  from  Figure  II- 1 that  the  stations  KIP  and  CTA 
are  almost  on  the  same  event- station  great  circle  path.  Referring  to  Table 
II- 3a,  the  azimuthal  angles  at  the  source  for  these  two  stations  differ  by  only 
1 degree,  and  there  are  four  events  (NTS/103/76,  NTS/314/76,  NTl/317/76, 
and  NT2/317/76)  which  have  Rayleigh  wave  recordings  at  these  two  stations. 
Therefore,  in  addition  to  obtaining  the  average  Rayleigh  wave  attenuation 
coefficients  by  the  isotropic- source  method  for  the  first  group  of  travel  paths, 
the  two- station  method  is  used  to  obtain  the  Rayleigh  wave  attenuation  coef- 
ficients particularly  for  the  travel  path  between  stations  KIP  and  CTA. 
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Figure  IV- 1 shows  the  Rayleigh  wave  amplitude  attenuation 
coefficients  from  10  to  50  seconds  calculated  at  5 second  increments,  and  ob- 
tained for  the  first  group  of  travel  paths.  In  this  figure,  the  average  Rayleigh 
wave  amplitude  attenuation  coefficients  k are  obtained  by  the  isotropic- 
source  method  for  the  travel  paths  between  the  NTS  and  the  stations  CTA, 

KIP,  and  MAT.  These  are  indicated  by  the  solid  dot  (•)  with  the  vertical  bar 

showing  the  standard  error  ±S  of  the  estimated  k value.  The  Rayleigh 

A 

wave  amplitude  attenuation  coefficients  obtained  by  the  two- station  method  for 
the  travel  path  between  the  stations  KIP  and  CTA  are  indicated  by  the  cross 
(x).  For  comparison  purposes,  the  average  Rayleigh  wave  attenuation  coef- 
ficients obtained  by  Mitchell  et  al.  (1976)  for  the  travel  paths  within  the  Pacific 
Ocean  are  indicated  by  the  solid  square  (■).  The  values  of  Mitchell  et  al. 

(1976)  agree  quite  well  with  our  values,  except  that  their  value  at  1 5 seconds 
is  significantly  higher.  Nevertheless,  all  their  values  lie  well  within  the 
standard  error  of  our  values.  Also,  the  agreement  between  the  values  for  the 
travel  path  from  the  station  KIP  to  the  station  CTA  and  the  average  values  for 
the  travel  paths  from  the  NTS  to  the  stations  CTA,  KIP,  and  MAT  is  good. 
Compared  to  the  values  of  Mitchell  et  al.  (1976),  our  values  exhibit  larger 
scatter  at  the  short  periods.  This  probably  is  due  to  our  using  fewer  obser- 
vation stations  (3  stations  versus  20  stations).  Also,  some  of  the  scatter  is 
due  to  variations  of  the  observed  station  spectral  amplitudes  from  anisotropic 
surface  wave  radiation  caused  by  not  accounting  for  tectonic  strain  releases 

of  the  NTS  events  in  using  the  isotropic  source  method.  The  above  facts  are 

- -4-1 

also  thought  to  be  responsible  for  the  larger  standard  error  (Ssl.SxlO  km  ) 

as  compared  to  the  variation  Idk^l  of  the  calculated  k.  value  by  the  two- 

^-4  -1 

station  method,  given  by  equation  (IV-4)  5:1.0  x 10  km  for 

D^-Dj  = 7610  km  between  KIP  and  CTA). 

Figure  IV- 2 shows  the  average  Rayleigh  wave  amplitude  atten- 
uation coefficients  from  10  to  50  seconds  at  5 second  increments.  These 
were  obtained  from  the  second  group  of  travel  paths.  These  values  are 
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Amplitude  Attenuation  Coefficient 


Travel  paths  from  NTS  to  ALPA,  NORS,  LASA, 
CHG,  TLO,  KON,  OGD,  ALQ,  ZLP,  and  MAI 


©Travel  paths  in  north  America  (Herrmann  & Mitchell,  1 975 ) 


5.  0 
4,  0 

^-1 

£ 3.  0 

'o  2.  0 
X 

< 1.0 

U 

0.  0 
-1.  0 

0 10  20  30  40  50 

Period  (sec) 


FIGURE  IV- 2 

AVERAGE  RAYLEIGH  WAVE  AMPLITUDE  ATTENUATION  FOR 
TRAVEL  PATHS  FROM  THE  NTS  TO  OBSERVATION  STATIONS 
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indicated  by  the  solid  dot  (•)  with  the  vertical  bar  showing  the  standard  error. 

These  values  indicate  that  the  average  Rayleigh  wave  amplitude  attenuation 

-4  -1 

for  the  second  group  of  travel  paths  decreases  rapidly  from  4.  1 x 10  km 

-4  - 1 

at  10  seconds  to  1, 8 x 10  km  at  20  seconds  and  stays  at  a lower  value 
-4  - 1 

(about  1.6x10  km  ) between  25  seconds  and  50  seconds.  For  comparison 
purposes,  this  figure  shows  the  average  Rayleigh  wave  attenuation  coefficients 
obtained  by  Herrmann  and  Mitchell  (1975)  for  the  travel  paths  within  North 
America,  and  the  Rayleigh  wave  attenuation  curve  given  by  Tryggvason  (1965) 
for  travel  paths  from  Novaya  Zemlya  to  stations  in  Europe,  Asia,  and  North 
America.  Herrmann  and  Mitchell's  values  are  shown  as  open  circles  (o) 
while  Tryggvason's  values  are  given  by  the  solid  curve.  Results  from  the 
comparison  among  these  three  sets  of  Rayleigh  wave  attenuation  values  can 
be  summarized  as  follows: 

• The  variations  of  the  k^  values  with  respect  to  period  are 
very  similar  among  these  three  sets.  They  all  indicate  that 
the  attenuation  is  highest  at  10  seconds,  decreases  rapidly  to 
a low  value  at  about  20  to  25  seconds,  and  maintains  a low 
value  with  very  little  change  between  25  and  50  seconds. 

• In  general,  among  these  three  sets  of  attenuation  values,  our 
values  are  the  highest  and  Tryggvason's  are  the  lowest.  How- 
ever, this  difference  is  thought  not  to  be  significant,  because: 
(1)  all  Herrmann  and  Mitchell's  values  fall  well  within  the 
standard  error  of  our  values  and  most  of  Tryggvason's  values 
also  lie  within  our  standard  error,  and  (2)  the  standard  er- 
rors of  these  three  sets  of  k . values  overlap, 

A 

• In  the  period  range  from  1 0 to  50  seconds,  the  standard  error 

-4  - 1 

of  our  values  is  about  1.  35  x 10  km  and  those  of  Tryg- 

-4 

gvason's,  and  Herrmann  and  Mitchell's  are  about  0.  8 x 10  km 
-4  -1 

and  2.1  x 10  km  , respectively.  The  smaller  standard 
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error  of  Tryggvason's  values,  as  compared  to  the  standard 
error  of  our,  can  probably  be  attributed  to  the  fact  that  (1) 
Tryggvason  used  a much  larger  number  of  stations  (14  sta- 
tions versus  our  10  stations),  and  (2)  the  travel  paths  involved 
in  Tryggvason's  estimation  are  more  uniform.  However,  the 
larger  standard  error  of  Herrmann  and  Mitchell's  values,  as 
compared  to  the  standard  error  of  ours,  seems  unreasonable 
since  Herrmann  and  Mitchell  used  a much  larger  number  of 
stations  than  we  did  in  the  estimation  of  values  (about  35 

stations  versus  our  10  stations),  and  their  method  did  account 
for  the  spectral  amplitude  variation  among  various  stations 
due  to  anisotropic  surface  wave  radiation.  Several  possible 
explanations  for  their  larger  standard  error  are:  (1)  the 
source  parameters  which  they  assumed  for  their  events  might 
not  be  correct,  while  accurate  prior  knowledge  of  the  event 
source  parameters  is  quite  critical  to  their  method,  and  (2) 
the  quality  of  their  data  might  not  be  uniformly  good  for  var- 
ious observation  stations, 

2.  Travel  Paths  from  EKZ  to  the  Observation  Stations 

Referring  to  Figure  III-2,  it  can  be  observed  that  almost  all  of 
the  travel  paths  encountered  from  EKZ  to  the  observation  stations  are  con- 
tinental or  mixed.  Referring  to  Table  II- lb  for  the  observation  stations 
available  to  the  selected  EKZ  events,  the  travel  paths  which  are  included  in 
the  estimation  of  the  Rayleigh  wave  amplitude  attenuation  are  those  from 
EKZ  to  stations  ALPA,  NORS,  LASA,  TLO,  EIL,  KON,  KIP,  ALQ,  ZLP, 
MAT,  GUM,  and  MAI.  The  travel  path  from  EKZ  to  station  CHG  is  intention- 
ally excluded  since  surface  waves  traveling  along  this  path  exhibit  an  abnormal 
group  delay  (as  discussed  in  Subsection  III-A)  and  abnormal  absorption  (as 
observed  by  Turnbull,  et  al,  , 1973).  The  isotropic-source  method  is  used  to 
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estimate  the  average  Rayleigh  wave  amplitude  attenuation  coefficients  from 
10  to  50  seconds  at  5 second  increments.  It  is  noticed  from  Figure  III- 2 that 
stations  KON  and  ZLP  almost  have  the  same  event- station  great  circle  path 
from  EKZ.  Referring  to  Table  II- 3b,  the  azimuthal  angles  at  the  source  for 
these  two  stations  differ  by  2.  1 degrees.  Therefore,  the  two- station  method 
is  applied  to  event  EKZ/704/76,  which  had  Rayleigh  wave  recordings  at  these 
two  stations,  to  calculate  the  Rayleigh  wave  attenuation  coefficients  along  the 
travel  path  between  stations  KON  and  ZLP.  This  travel  path  is  mostly  oce- 
anic, across  the  Atlantic  Ocean. 

Figure  IV- 3 presents  the  obtained  Rayleigh  wave  amplitude  at- 
tenuation coefficients  for  the  travel  path  between  stations  KON  and  ZLP.  They 
are  shown  by  solid  dots  (•).  Also  included  in  this  figure  are  the  Rayleigh  wave 
amplitude  attenuation  coefficients  previously  obtained  for  the  travel  path  be- 
tween stations  TLO  and  ZLP  (Sun,  1976)  and  those  obtained  by  Tsai  and  Aki 
(1969)  for  the  travel  path  between  station  BEC  in  Bermuda  and  station  SDB  in 
South  Africa.  They  are  shown  by  crosses  (x)  and  by  solid  squares  (■),  res- 
pectively. These  three  travel  paths  are  mostly  oceanic,  across  various  parts 
of  the  Atlantic  Ocean,  and  their  Rayleigh  wave  attenuation  coefficients  shown 
in  Figure  IV- 3 are  all  obtained  by  the  two- station  method.  The  variations 
among  these  three  sets  of  attenuation  coefficients  are  small,  except  that  Tsai 
and  Aki's  (1969)  values  show  larger  scattering  and  a higher  attenuation  at  1 5 
seconds. 

Figure  IV-4  shows  the  average  Rayleigh  wave  amplitude  at- 
tenuation coefficients  for  the  travel  paths  from  EKZ  to  the  observation  sta- 
tions. These  values  are  indicated  by  the  solid  dots  (•)  with  vertical  bars 
showing  the  standard  error.  These  values  indicate  that  the  average  Rayleigh 

wave  amplitude  attenuation  for  these  travel  paths  decreases  rapidly  from 

-4  -1  -4  -1 

4.  2 X 10  km  at  10  seconds  to  1. 6 x 10  km  at  20  seconds  and  decreases 

gradually  with  small  variations  from  25  to  50  seconds.  In  this  figure,  again 
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• Travel  path  between  KON  and  ZLP 


FIGURE  IV- 3 

RAYLEIGH  WAVE  AMPLITUDE  ATTENUATION  FOR 
TRAVEL  PATHS  IN  THE  ATLANTIC  OCEAN 
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1 Travel  paths  from  EKZ  to  ALFA,  NORS,  LASA,  TLO, 

' EIL,  KON,  KIP,  ALQ,  ZLP,  MAT,  GUM,  and  MAI 

o Travel  paths  in  north  America  (Herrmann  & Mitchell,  19(5) 


Travel  paths  in  Europe,  Asia,  and  north  America 
(Tryggvason,  1965) 
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for  comparison  purposes,  the  average  Rayleigh  wave  amplitude  attenuation 
coefficients  obtained  by  Herrmann  and  Mitchell  (1975)  and  those  given  by 
Tryggvason  (1965)  are  shown  by  open  circles  (o)  and  solid  curves,  respective- 
ly. Results  from  the  comparison  among  these  three  sets  of  Rayleigh  wave 
amplitude  attenuation  values  can  be  summarized  as  follows: 

• The  variations  of  the  values  with  respect  to  the  period  are 

quite  similar  among  these  three  sets.  They  all  indicate  that 
the  attenuation  is  the  highest  at  1 0 seconds,  decreases  rapidly 
toward  20  seconds  and  stays  at  a lower  value  with  small  varia- 
tions from  25  to  50  seconds. 

• Herrmann  and  Mitchell's  values  agree  very  well  with  our  values 
and  all  their  values  lie  well  within  the  standard  error  of  our 
values.  Tryggvason's  values  are  somewhat  lower  than  ours, 
yet  most  of  Tryggvason's  values  fall  within  our  standard  error. 
Thus,  this  difference  in  the  attenuation  values  is  not  thought  to 
be  too  significant. 

• For  the  period  range  from  10  to  50  seconds,  the  standard  error 

-4  - 1 

of  our  value  is  about  1.0x10  km  . The  remarks  made 

A 

in  Subsection  IV-B-1,  regarding  the  standard  errors  of  the 
three  methods,  also  apply  here. 
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SECTION  V 


SOURCE  PARAMETER  ESTIMATES  OF  THE  SELECTED 
NTS,  EKZ,  AND  PNE  EVENTS 


Many  underground  nuclear  explosions  generate  horizontally 
polarized  seismic  shear  waves  (SH  and  Love)  along  with  P,  SV,  and  Ray- 
leigh waves.  Theoretically,  an  explosive  source  in  a horizontally  layered, 
homogeneous,  isotropic  medium  should  not  generate  SH  or  Love  waves.  The 
significance  of  the  magnitude  of  the  SH  and  Love  waves  radiated  from  explo- 
sions has  been  demonstrated  by  many  investigators  (Press  and  Archambeau, 
1962;  Brune  and  Pomeroy,  1963;  Toksbz,et  al.  , 1965;  and  Lambert,  et  al.  , 
1972).  The  possible  mechanisms  which  might  be  responsible  for  these  SH  and 
Love  waves  have  also  been  discussed  by  them.  These  mechanisms  are  wave 
type  conversion,  near- source  irregularity,  and  tectonic  strain  release.  The 
first  two  mechanisms  have  been  ruled  out  by  simple  comparison  of  various 
observations  from  different  explosions  in  the  same  area  (Prune  and  Pomeroy, 
1963;  and  Toksbz.et  al.  , 1965).  The  tectonic  strain  release  associated  with 
the  explosion  is  thought  to  be  the  most  plausible  mechanism  for  the  genera- 
tion of  the  SH  and  Love  waves  by  the  explosion  event  (Toksoz.et  al,  , 1965), 

A nuclear  explosion  is  not  an  ideal  point  source.  The  shock 
wave  generated  vaporizes,  melts,  crushes,  and  cracks  the  rocks  out  to  some 
distance,  depending  on  the  yield  of  the  explosion  and  the  strength  of  the 
medium.  The  cracks  caused  by  the  explosion  would  not  generate  any  SH  or 
Love  waves  unless  there  were  differential  movement  along  them.  In  a medium 
which  is  not  pre- stressed  the  cracks  and  the  movement  along  cracks  would  be 
distributed  uniformly,  with  radial  symmetry  around  the  point  of  explosion. 
Therefore,  one  would  not  expect  an  SH  or  Love  wave  radiation  pattern  which 
could  be  represented  by  a dipole  or  a double-couple  source.  In  a pre- stressed 
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medium,  an  explosion  may  release  a part  of  the  existing  strain  energy  by  one 
or  more  of  the  following  three  mechanisms.  First,  the  intense  stress  waves 
near  the  explosion  may  trigger  an  earthquake;  second,  the  strain  energy  may 
be  released  by  directional  cracking;  and  third,  the  strain  energy  may  be  re- 
leased by  the  creation  of  a cavity. 

In  this  section,  our  goal  is  not  to  verify  the  existence  of  the 
SH  and  Love  waves  from  the  explosion  nor  to  demonstrate  the  possible  mech- 
anisms which  generate  these  shear  waves.  Rather,  in  this  section,  our  spec- 
tral fitting  technique  is  designed  to  take  into  account  any  possible  tectonic 
strain  release  associated  with  an  explosion.  The  technique  is  applied  to 
selected  NTS,  EKZ,  and  PNE  explosion  events  (as  listed  in  Table  II- 1). 

Source  parameters  of  these  events  will  be  determined  from  the  observed  far- 
field  surface  wave  spectra.  Hopefully,  source  parameter  estimates  of  these 
events  will  lead  to  a better  understanding  of  the  physics  of  the  explosions  from 
NTS  and  EKZ. 

If  there  is  tectonic  strain  release,  P,  SV,  SH , Love,  and  Ray- 
leigh waves  are  radiated.  These  will  be  added  to  the  explosion -gene rated  P, 

SV,  and  Rayleigh  waves.  Although  the  SH  and  Love  waves  will  be  due  to  the 
tectonic  strain  release  only,  the  P,  SV,  and  Rayleigh  waves  will  contain  frac- 
tions due  to  both  sources.  It  is  not  possible  to  separate  the  fraction  of  energy 
due  to  the  explosion  alone  unless  the  space  and  time  functions  of  both  sources 
are  known  exactly.  Since  this  information  is  not  available,  we  resort  to  some 
indirect  techniques  and  make  certain  assumptions,  which  will  be  discussed 
later  on. 

In  the  following,  a procedure  of  source  parameter  estimation 
using  far-field  surface  wave  data  will  be  described  briefly  in  Subsection  A. 
Subsection  B will  present  the  layered  half- space  earth  models  used  to  analyze 
the  NTS,  EKZ,  and  PNE  events.  In  Subsection  C,  the  source  parameter  esti- 
mates of  the  selected  NTS,  EKZ,  and  PNE  events  will  be  presented  and  the 
results  from  the  examination  of  these  source  parameter  estimates  will  be  given. 
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A. 


A PROCEDURE  OF  SOURCE  PARAMETER  ESTIMATION  USING  FAR- 
FIELD  SURFACE  WAVE  DATA 


A procedure  to  obtain  the  seismic  source  parameter  estimation 
using  far-field  surface  wave  data  is  schematically  shown  in  Figure  V-).  The 
field  tapes  from  observation  stations  are  edited  to  obtain  the  surface  wave 
time  traces  for  a particular  event.  The  edited  time  traces  are  analyzed  using 
the  ISPS  to  obtain  the  demultipathed  surface  wave  spectra,  as  discussed  by 
Sun  and  Shaub  (1977).  The  observed  spectra  are  corrected  for  instrument 
response,  geometric  spreading,  and  travel  path  attenuation.  Finally,  the 
corrected  spectra  are  input  to  an  amplitude  spectral  fitting  process  to  obtain 
estimates  of  the  source  parameters.  The  first  three  blocks  in  Figure  V-1 
are  pertinent  to  the  data  handling  and  have  been  discussed  in  Sections  II  and 
III,  In  the  following,  the  amplitude  spectral  fitting  process  will  be  briefly 
described. 

The  amplitude  spectra  fitting  process  employs  the  modified 
Tsai's  method  for  the  estimation  of  seismic  source  parameters  (Tsai  and 
Aki,  1970;  Tsai  and  Shen,  197Z:  Turnbull, et  al.  , 1973).  The  modified  Tsai's 
method  is  an  exhaustive  searching  scheme  to  find  the  best  fit  between  the  ob- 
served and  the  theoretically  calculated  surface  wave  amplitude  spectra.  Or- 
iginally, the  modified  Tsai's  method  was  developed  to  obtain  a source  param- 
eter estimate  of  an  earthquake  modeled  by  a double-couple  source.  The 
model  yields  five  source  parameters;  seismic  moment  (M),  focal  depth  (h), 
dip  (<J),  slip  (X).  ?nd  strike  (0)  angle.  For  an  explosion  event  with  associated 
tectonic  strain  release,  the  theoretical  source  model  which  is  used  here  com- 
bines a point  explosive  source  and  a point  double-couple  source.  Theoretical 
surface  wave  spectra  are  generated  by  this  combined  source  in  a layered  half- 
space earth  (Harkrider,  1964;  Ben-Menahem  and  Harkrider,  1964;  Toksdz, 
et  al.  , 1965).  A point  explosive  source  models  an  ideal  explosion  and  radiates 
only  isotropic  Rayleigh  waves,  while  a point  double-couple  source  models  the 
explosion-associated  tectonic  strain  resease  and  radiates  anisotropic  Rayleigh 
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and  Love  waves.  Therefore,  the  resultant  Rayleigh  wave  radiated 
from  the  combined  source  is  the  vector  sum  of  the  Rayleigh  wave  from  the 
point  explosion  source  and  that  from  the  point  double-couple  source.  The  re- 
sultant Love  wave  from  the  combined  source  comes  entirely  from  the  point 
double-couple  source.  The  resultant  amplitude  spectra  can  be  expressed  as 
follows : 


M^LR^ 
E E 


(h.<J 


+ M LR  (h)  .. 

A.  A IJ 


where 


LR_  + F LR^,  . . 
E X ij 


LQ„(h,6,X,0)  .. 
E ij 


(V-1) 


Y^.  and 
iJ 


LR, 


and  LQ_ 
E 


LR 


X 


M and  M 
£ X 

F 


are  the  theoretical  Rayleigh  and  Love  wave 
amplitude  spectra  of  the  combined  source, 
respectively,  for  the  i-th  period  at  the  j-th 
station, 

are  the  theoretical  Rayleigh  and  Love  wave 
spectra  of  the  double-couple  source,  res- 
pectively, for  a seismic  moment  of  one  dyne- 
cm, 

is  the  theoretical  Rayleigh  wave  spectrum  of 
the  point  explosive  source  for  a seismic 
moment  of  one  dyne-cm, 

are  the  moments  of  the  double-couple  source 
and  the  explosive  source,  respectively,  and 
is  the  relative  strength  of  the  explosive  source 
with  respect  to  the  double-couple  source,  i.  e.  , 
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It  is  assumed  here  that  the  explosion-associated  tectonic  strain  release  takes 
place  in  the  vicinity  of  the  explosion,  i.  e.  , the  depths  of  the  explosive  source 
and  the  double-couple  source  are  the  same.  Therefore,  for  an  explosion  with 
tectonic  strain  release,  there  are  six  parameters  involved  in  the  theoretical 
source  model.  These  are  focal  depth  (h),  dip  angle  (i),  slip  angle  (X),  strike 
angle  (0),  seismic  moment  of  the  double-couple  source  (M„),  and  the  relative 
strength  of  the  explosion  (F).  The  geometry  of  the  double-couple  source  which 
models  the  explosion-associated  tectonic  strain  release  is  shown  in  Figure  V-Z 

Let  s ((5,\,0,h,  F)  be  a member  in  a predefined  set  of  source 
P 

parameters  S((5 , X , 0 , h,  F).  The  exhaustive  search  in  the  amplitude  spectral 
fitting  process  tests  each  member  in  S and  obtains  the  member  which  yields 
the  best  fit  between  the  observed  and  the  theoretically  calculated  surface  wave 
amplitude  spectra.  The  entire  set  of  the  source  parameters  S can  be  con- 
veniently defined  by  the  upper  and  lower  bounds  of  each  source  parameter  with 
discrete  increment;  that  is 


h < h<  h 
1 m 

6 < 6 

1 k 

Aj<  X < 

4>  < <l>  < d> 

1 ■ n 

F < F < F 

1 q 


with  increment  of  dh 
with  increment  of  46 
with  increment  of  4 A 
with  increment  of  40 
with  increment  of  4F 


where  m,  k,  I,  n,  and  q are  positive  integers.  Thus,  the  entire  set  S is  a 
discrete  set  consisting  of  m*  k*  f • n*q  members.  In  the  amplitude  spectral 
fitting  process  which  is  applied  here,  the  following  values  of  the  source 
parameters  are  used: 


h = 0,5,  1,5,  2.5,  4.0,  6.0,  8.0,  10,0  km; 
30°<  6 i 9o”,  46  = 10°:  -90°<  A < 90°,  4A  = 30°; 
0°S  0 < 175°,  40  = 5°;  0.  25  < F<  3.  0,  4F  = 0,  25. 
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FIGURE  V-2 

GEOMETRY  OF  THE  DOUBLE- COU RLE  SOURCE  WHICH  MODEI 
THE  EXPLOSION-ASSOCIATED  TECTONIC  STRAIN  REl.EAM 
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The  non-uniform  increment  in  h is  dictated  by  the  individual  layer  thickness 
of  the  layered  half- space  earth  model  which  is  chosen  for  the  source  region. 

R L 

Let  X..  and  X..  denote  the  observed  Rayleigh  and  Love  wave 
ij  iJ 

spectral  amplitudes,  respectively,  for  the  i-th  period  at  the  j-th  station.  The 
difference  between  the  observed  spectra  and  the  corresponding  theoretical 
spectra,  which  will  be  called  the  residual  €,  is  defined  as: 


£(Sp)  = 


??K- 


y!^)^  + (X^  - Y 


J • 


Substituting  (V-l)  into  (V-2),  we  have 


(V-2) 


€(s  ) 
P 


LRe  + flr^ 


+ (X^  - M 


E 


(V-2a) 


In  the  amplitude  spectra  fitting  process,  first  the  moment  is  estimated 

by  minimizing  the  residual  €(s  ) for  all  i and  j , or  explicitly: 


EE 

j i 

X^ 

LRe+FLRx 

. . + X.^ 
ij 

EE 

j i 

LR 

2 

. . + 
ij 

'■“e  I ij] 

In  other  words,  for  a particular  member  Sp  in  the  set  of  source  parameters 
S,  the  optimum  moment  is  first  found  from  equation  (V-3).  Then,  the  result- 
ant ME(Sp)  is  used  in  equation  (V-2a)  to  compute  the  minimum  residual  ^(Sp) 

for  this  particular  member  s . The  minimum  residual  € (s  ) is  computed 

P P 

for  each  of  the  members  in  S . This  € value  is  used  as  a convenient  cri- 
terion for  judging  the  goodness  of  the  fit  between  the  observed  and  the  theo- 
retically calculated  spectra  for  each  member  s . In  the  end,  the  source 


parameter  estimation  of  a given  event  will  be  the  member  s 


optimal 


whose 


theoretical  surface  wave  spectra  fits  most  closely  the  observed  spectra;  that 


is: 
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^all  optimal  ^ Minimum  ft {.  )1  for  all 
min  p L P J 


8 in  S. 
P 


(V-4) 


In  order  to  attain  a greater  degree  of  confidence  in  using  the 
source  parameter  estimation  based  on  the  minimum  residual  criterion,  a de- 
tailed knowledge  of  the  residual  distribution  about  the  absolute  minimum  given 
by  equation  (V-4)  would  be  helpful  to  guard  against  a possible  spurious  mini- 
mum generated  by  poor  quality  data.  Therefore,  in  the  amplitude  spectral 
fitting  process  used  here,  the  distributions  of  f (8p)  with  respect  to  the 
source  parameters  are  generated  as  follows. 

In  the  fitting  process,  the  focal  depth  is  varied  last,  because 
the  focal  depth  (h)  will  affect  the  shape  of  the  theoretical  surface  wave  ampli- 
tude spectra  the  most,  especially  for  Rayleigh  waves  (Turnbull,et  al.  , 1974). 

There  is  one  c^i.  at  one  focal  depth  h.  , that  is 
min  ^ 1 


min 


= Minimum  , X , 0 , h^  F)J| 


for  all  d,X  ,<t>  and 
F at  depth  h^  . 


(V-5) 


The  variation  of  6 i.  with  respect  to  the  focal  depth  is  obtained  and  a stand- 


ard value  f , say 


[■«'“  1 
L minJ 


^ , ,,  is  determined  from  those  f for 

standard  min 


the  purpose  of  picking  appropriate  members  in  the  set  of  source  parameters 
S to  form  a sample  set  . In  the  entire  set  of  source  parameters  S , those 
members  whose  €(8^)'s  are  less  than  €^i.  , , , are  included  in  the 

sample  set  S^.  These  form  the  distributions  of  €{Sp)  with  respect  to  the 
source  parameters.  These  residual  distributions  for  the  various  source 
parameters  will  be  more  easily  understood  from  the  results  which  will  be 
presented  later  on  lor  the  source  parameter  estimations  of  the  selected  NTS, 


EKZ,  and  PNE  events. 


B.  EARTH  MODELS 

In  order  to  generate  the  theoretical  surface  wave  spectra  which 
are  needed  in  the  amplitude  spectral  fitting  process,  a layered  half-space, 
which  models  the  source  region  earth  structure,  has  to  be  obtained  first. 
Three  earth  models  are  used  for  the  selected  NTS,  EKZ,  and  PNE  events, 
respectively. 

For  the  NTS  events,  two  structures  appear  to  be  appropriate 
for  the  source  mechanism  study.  One  is  the  35CM2  earth  model  of  Alexander 
(1963)  for  the  events  occurring  near  (37°N,  116°W),  and  the  other  is  the  mod- 
ified Eaton's  model  (Toksbz,et  al.  , 1965)  for  the  events  near  (39°N,  118°W). 
Since  our  selected  NTS  events  all  took  place  near  (37°N,  116°W),  the  layered 
half- space  earth  model  for  our  selected  NTS  events  was  constructed  from  the 
35CM2  model.  Basically,  the  35CM2  model  is  a simple  structure  consisting 
of  about  2-3  km  of  'sedimentary'  material  overlaying  a rather  uniform  section 
of  'granite'  material  about  23  km  thick,  which  in  turn  overlies  an  intermediate 
section  about  25  km  thick.  A maximum  propagation  velocity  occurs  at  a depth 
of  about  50  km  with  a gradual  decrease  in  velocity  at  greater  depth  corres- 
ponding to  the  beginning  of  the  Gutenberg  low-velocity  zone  centered  at  a depth 
of  around  125-150  km.  To  estimate  source  parameters  from  long-period  sur- 
face waves,  we  need  the  structure  profile  of  the  earth  model  down  to  a depth 
of  around  400-500  km.  Since  the  35CM2  model  has  a similar  structure  pro- 
file as  the  Gutenberg  model  for  any  depth  greater  than  about  50  km,  our  NTS 
earth  model  merges  the  35CM2  model  with  the  Gutenberg  model  at  a depth  of 
150  km.  The  structure  profile  of  our  NTS  earth  model  is  given  in  Table  V-1, 
and  theoretical  dispersions  for  the  Rayleigh  and  the  Love  wave  are  given  in 
Table  V-2. 

For  the  selected  EKZ  and  PNE  events,  the  earth  models  are 
constructed  from  the  information  given  by  Rodriguez  (1969).  The  locations 
of  our  selected  EKZ  and  PNE  events  are  shown  on  Figure  V-3,  in  relation  to 
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TABLE  V-1 

STRUCTURE  PROFILE  OF  THE  NTS  EARTH  MODEL  (37°N,  116°W) 


Thickness  (km) 

Vp  (km/ sec ) 

Vg  (km/sec) 

p(g/cm^) 

2 

3.  80 

1. 75 

2.  20 

13 

6,  10 

3.  60 

2.  80 

10 

6.  20 

3.  70 

2.  80 

11 

7.62 

4.  07 

3.  20 

14 

7.  62 

4.  07 

3.  25 

16 

8.  10 

4.  62 

3.  32 

20 

8.  00 

4.  53 

3.  40 

20 

7.94 

4.  46 

3.40 

45 

7.  80 

4.  40 

3.  40 

25 

7.  89 

4.  35 

3.43 

25 

7.  98 

4.  38 

3.  46 

25 

8.  10 

4.42 

3.48 

25 

8.  21 

4.46 

3.  50 

50 

8.  38 

4.  54 

3.  53 

50 

8,  62 

4.  68 

3.  58 

50 

8,  87 

4,  85 

3.  62 

50 

9.  15 

5.  04 

3.  69 

50 

9.  45 

5.  21 

3.  82 

50 

9.  88 

5.45 

4.  01 

Rayleigh  Wave  (LR)  Love  Wave  (LQ) 


TABLE  V-2 


P 


SURFACE  WAVE  DISPERSION  FOR  THE  NTS  EARTH  MODEL 

(37°N,  llb°W) 


T (sec) 

C (km/ sec ) 

V (km/sec) 

k (km  ) 

60 

4.  351 

4.  107 

0. 0241 

55 

4.  328 

4.  069 

0. 0264 

50 

4.  300 

4.  018 

0. 0292 

a 

45 

4.  266 

3.  948 

0. 0327 

4) 

40 

4.  222 

3.  856 

0. 0372 

35 

4.  165 

3.  747 

0. 0431 

> 

30 

4.  089 

3.  632 

0,  0512 

0 

25 

3.  993 

3.  527 

0. 0629 

20 

3.  878 

3.  436 

0. 0810 

15 

3.  742 

3.  342 

0.  112 

10 

3.  564 

3.  161 

0.  176 

60 

3.  981 

3.  877 

0. 0263 

55 

3.  971 

3.  855 

0.  0288 

50 

3.  959 

3.  817 

0.  0318 

45 

3.  941 

3.  755 

0.  0354 

V 

> 

Id 

40 

3.  914 

3,  661 

0.  0401 

35 

3.  871 

3.  526 

0.  0464 

OD 

4) 

30 

3.  802 

3.  353 

0.  0551 

«j 

25 

3.  698 

3.  170 

0.  0680 

cc 

20 

3.  557 

3.  020 

0.  0883 

15 

3.  386 

2.  916 

0.  124 

10 

3.  203 

2.  886 

0.  196 

Tj  Period,  C:  Phase  Velocity,  V:  Group  Velocity,  k;  Wavenumber 
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Latitude 


LOCATIONS  OF  THE  SELECTED  EKZ  AND  PNE  EVENTS,  AND  LOCAL  GEOLOGY 


major  geological  features.  The  selected  EKZ  events  are  in  the  Balkhash- 
Chingiz  Foldbelt,  while  the  selected  PNE  events  are  in  the  north  Caspian 
Basin.  For  EKZ,  the  earth  model  is  constructed  from  a representative  sub- 
surface crustal  structure  of  the  Balkhash- Chingiz  Foldbelt;  and,  for  PNE, 
the  earth  model  is  constructed  from  a representative  subsurface  crustal 
structure  of  the  north  Caspian  Basin.  The  crustal  cross  sections  along  the 
broken  lines  XIV-XIV  and  C-C  are  shown  in  Figures  V-4  and  V-5,  res- 
pectively, giving  the  major  structure  parameters.  Since  these  two  crustal 
cross  sections  only  go  down  to  the  depth  around  50  km,  we  have  again  con- 
structed our  EKZ  and  PNE  earth  model  by  merging  the  crustal  cross  sections, 
given  in  Figures  V-4  and  V-5,  respectively,  with  the  Gutenberg  model  at  the 
depth  of  50  km.  The  structure  profile  of  our  EKZ  earth  model  is  given  in 
Table  V-3,  and  the  theoretical  dispersions  for  the  Rayleigh  and  Love  waves 
are  given  in  Table  V-4.  Similar  information  for  our  PNE  earth  model  is 
given  in  Tables  V-5  and  V-6. 

C.  SOURCE  PARAMETER  ESTIMATIONS 

The  observed  surface  wave  spectra  of  the  selected  NTS,  EKZ, 
and  PNE  events  have  been  obtained  in  the  manner  described  in  Section  II.  Us- 
ing the  appropriate  layered  half- space  earth  model,  the  amplitude  spectral 
fitting  process,  as  discussed  in  Subsection  V-A,  is  applied  to  the  observed 
surface  wave  spectra  of  these  selected  events  to  obtain  source  parameter 
estimates.  In  the  following,  the  results  of  the  source  parameter  estimations 
will  be  presented  separately  for  selected  NTS,  EKZ,  and  PNE  events,  and 
comparisons  of  these  results  will  be  made. 

1.  The  NTS  Events 

For  NTS  events,  the  estimations  of  the  source  parameters 
based  on  the  minimum- residual  criterion,  as  defined  in  equation  (V-4),  are 
listed  in  Table  V-7.  Except  as  noted  in  the  table,  the  source  parameter 
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TABLE  V-4 


SURFACE  WAVE  DISPERSION  FOR  THE  EKZ  EARTH  MODEL 
(BELKHASH-CHINGIZ  FOLDBELT) 


T (sec) 

C (km/sec) 

V (km /sec) 

k (km  ^ ) 

60 

4.  297 

4.  010 

0. 0244 

55 

4.  270 

3.  960 

0, 0268 

50 

4.  236 

3,  892 

0, 0297 

a 

45 

4.  193 

3.  801 

0. 0333 

0) 

> 

40 

4.  138 

3.  687 

0. 0380 

35 

4.  065 

3.  558 

0.  0442 

> 

o 

30 

3.  972 

3.  434 

0. 0527 

►J 

25 

3,  859 

3.  335 

0. 0651 

20 

3,  732 

3.  271 

0. 0842 

15 

3.  597 

3.  227 

0.  116 

10 

3.  455 

3.  182 

0.  182 

60 

3.  929 

3,  809 

0.  0267 

55 

3.  918 

3.  793 

0, 0292 

K 

50 

3.  905 

3.  762 

0. 0322 

45 

3.  888 

3.  708 

0. 0359 

> 

OJ 

40 

3.  862 

3.  617 

0.  0407 

X 

35 

3.  819 

3.  477 

0.  047 

00 

30 

3.  749 

3,  281 

0. 0559 

>s 

cx 

25 

3.  637 

3.  072 

0.  0691 

20 

3.  487 

2.  940 

0. 0901 

15 

3.  321 

2,  884 

0.  126 

■ 1 

10 

- - -- 

3.  147 

2.  816 

0.  200 

T:  Period,  C:  Phase  Velocity,  V;  Group  Velocity,  k:  Wavenumber 


V-18 


TABLE  V-5 


STRUCTURE  PROFILE  OF  THE  PNE  EARTH  MODEL 
(NORTH  CASPIAN  BASIN) 


Thickness  (km) 

Vp  (km/ sec) 

Vg  (km/sec) 

3 

P(g/cm  ) 

8 

6.  14 

3.  55 

2.  74 

16 

6,  75 

3.  90 

3.  05 

16 

7.  60 

4.  40 

3.  05 

10 

8,  20 

4.  65 

3.  32 

10 

8,  17 

4.  62 

3.  34 

10 

8.  14 

4.  57 

3.  35 

10 

8.  10 

4.  51 

3.  36 

10 

8.  07 

4,  46 

3.  37 

10 

8.  02 

4.  41 

3.  38 

25 

7.  93 

4.  37 

3.  39 

25 

7.  85 

4.  35 

3.  41 

25 

7.  98 

4.  38 

3.  43 

25 

8.  10 

4,  42 

3.  48 

25 

8,  21 

4.  46 

3,  50 

50 

8.  38 

4.  54 

3.  53 

50 

8.  62 

4.  68 

3.  58 

50 

8.  87 

4.  85 

3.  62 

50 

9.  15 

5.  04 

3.  69 

50 

9.  45 

5.  21 

3.  82 

50 

9.  88 

5.  45 

4.  01 

TABLE  V-6 


SURFACE  WAVE  DISPERSION  FOR  THE  PNE  EARTH  MODEL 
(NORTH  CASPIAN  BASIN) 


T (sec) 

C (km/sec) 

V (km/sec) 

k (krn'S 

60 

4.  434 

4.  291 

0. 0236 

55 

4.  421 

4.  286 

0. 0258 

50 

4.  408 

4.  277 

0. 0285 

o 

45 

4.  394 

4.  263 

0.  0318 

40 

4.  378 

4.  237 

0. 0359 

> 

rd 

35 

4,  357 

4.  191 

0.  0412 

30 

4.  326 

4.  106 

0.  048 

> 

0 

hJ 

25 

4,  275 

3.  967 

0. 0588 

20 

4.  189 

3.  806 

0, 0750 

15 

4,  065 

3.  683 

0.  103 

10 

3.  91^ 

3,  61 1 

0.  161 

60 

4.  022 

3.  990 

0.  026 

55 

4.  019 

4.  000 

0.  028 

50 

4,  018 

4.  005 

0.  0313 

CC 

45 

4,  017 

4,  000 

0. 0348 

0) 

> 

40 

4,  014 

3.  980 

0. 0391 

35 

4.  007 

3.  937 

0. 0448 

£ 

QD 

30 

3.  991 

3.  854 

0. 0525 

td 

25 

3.  955 

3.  713 

0. 0636 

a 

20 

3.  880 

3.  514 

0. 0810 

15 

3,  750 

3.  372 

0.  1 12 

10 

3.  581 

3.  267 

0.  175 

T:  Period,  C:  Phase  Velocity,  V;  Croup  Velocity,  k;  Wavenumber 
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ESTIMATIONS  OF  SOI  ROE  PARAMETERS  OBTAINED  BY  AMPLITUDE 
SPECTRAL  FITTING  BASED  ON  THE  MINIMUM- RESIDUAL 
CRITERION:  SELECTED  NTS  EVENTS 
(PAGE  2 OF  2) 


observation  spectra  at  the  SDCS  stations. 


estimates  of  these  events  are  obtained  without  including  the  observed  spectra 
at  the  SDCS  stations.  This  is  due  to  the  non-uniform  data  quality  observed 
at  the  SDCS  stations  as  discussed  in  Subsection  III-B.  The  depth  estimates 
of  these  events  are  all  very  shallow;  0.  5 km  except  2.  5 km  for  the  event 
NT2/204/76.  The  very  shallow  depths  of  these  events  are  to  be  expected, 
since  these  events  are  presumed  underground  explosions.  The  dip  and  slip 
angle  estimates  of  the  tectonic  strain  releases  associated  with  these  events 
can  be  generally  grouped  into  two  kinds:  nearly  strike- slip  for  10  events  and 
purely  dip- slip  for  6 events.  These  dip  and  slip  angle  estimates  are  quite 
different  from  the  assumption  made  by  other  investigators  (Toksoz.et  al.  , 
1965;  Lambert,et  al.  , 1972;  Mitchell,  1975)  that  the  tectonic  strain  release 
associated  with  the  NTS  explosion  has  the  vertical  strike- slip  type  of  mech- 
anism (6=  90°,  X = 0°).  Even  for  those  10  events  with  nearly  strike-slip  type 
of  mechanism  for  their  associated  tectonic  strain  releases,  there  are  varia- 
tions in  their  dip  and  slip  angle  estimates  (i  = 60°to  90°,  X=  -30°to30°).  How- 
ever, with  just  a few  exceptions,  the  strike  angle  estimates  are  quite  close 
to  the  N-S  direction  as  assumed  by  some  authors  (Toksbz,et  al.  , 1965;  Lam- 
bert,et  al.  , 1972).  The  estimates  of  the  relative  strength  of  the  explosive 
source  and  the  double-couple  source,  i.  e.  , the  F values  (M  /N'  .,),  vary 
from  1.  00  to  1.  75  among  these  events. 

Among  these  16  selected  NTS  events,  there  are  six  events 
which  can  be  grouped  into  three  pairs,  based  on  their  similar  locations  and 
times  of  detonations  (see  Table  II- 1).  These  three  pairs  are  NTl/603/75 
and  NT2/603/75,  NTl /31 7/76  and  NT2/31 7/76,  and  NTl /204/76  and 
NT2/204/76.  The  locations  of  two  events  in  each  pair  were  very  close  and 
their  detonations  were  about  20  to  30  minutes  apart.  The  source  param- 
eter estimates  of  tectonic  strain  releases  associated  with  two  events  in 

the  first  two  pairs  are  very  similar.  The  larger  event  (in  terms  of  its  m ) 

b 

of  each  pair  releases  a larger  amount  of  tectonic  strain  energy  (in  terms  of 
the  estimated  seismic  moment  of  the  tectonic  strain  release).  This  probably 
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indicates  that  the  two  events  in  each  of  the  first  two  pairs  had  similar  mech- 
anisms in  releasing  the  tectonic  strain  energy.  However,  the  source  param- 
eter estimates  of  the  tectonic  strain  releases  associated  with  the  two  events 
in  the  last  pair,  i.  e.  , NTl  /204/76  and  NT2/204/76,  indicate  quite  a different 
situation.  The  depth  estimate  of  the  event  NT2/204/76  is  2.  5 km  which  is 
quite  unusual  for  an  explosion  with  an  m^  of  5.  7.  The  dip  and  slip  estimates 
of  the  tectonic  strain  releases  associated  with  these  two  events  are  quite  dif- 
ferent and  the  event  NT2/204/76  has  a much  larger  seismic  moment  estimate 
of  the  tectonic  strain  release.  These  differences  in  source  parameter  esti- 
mates well  might  have  been  caused  by  a problem  with  the  observed  spectra, 
since  they  were  only  available  at  SDCS  stations  whose  observed  spectra  are 
questionable  as  discussed  in  Subsection  III-B.  Alternatively,  these  differences 
might  be  due  to  a real  physical  reason  if  NT2/204/76  released  the  tectonic 
strain  energy  by  actually  trigerring  an  earthquake  at  a depth  greater  than  the 
shot  depth. 

The  residual  distributions  with  respect  to  the  estimated  depth, 
dip  angle,  slip  angle,  and  strike  angle  are  shown  in  Figures  V-6  through 
Figure  V-21,  together  with  the  spectral  fit  for  each  event.  The  residual  dis- 
tribution for  one  particular  source  parameter,  say  depth,  consists  of  the 
variation  of  the  minimum  residual  values  with  respect  to  the  depths  tested  in 
the  spectral  fitting  process  and  the  percentage  of  occurrances  at  the  various 
tested  depths.  For  each  event,  the  residual  distributions  are  given  in  part 
(a)  of  each  figure.  The  best  spectral  fits  of  Rayleigh  and  Love  waves  (as  de- 
fined by  equation  (V-4))  are  shown  in  part  (b).  In  these  figures,  the  residual 

values  are  normalized  with  respect  to  and  are  indicated  as 

min  p 

the  relative  minimum  residuals  (R.  M.  R.  ).  The  observed  spectra  are  in- 
dicated by  crosses  (x)  and  the  best  fit  theoretical  spectra  are  represented 
by  open  circles  (o)  connected  by  straight  lines.  When  there  are  no  observed 
spectra  for  either  the  Rayleigh  or  Love  wave,  the  crosses  are  plotted  on  top 
of  the  corresponding  open  circles.  In  general,  the  agreement  between 
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FIGURE  V-6 

RESULTS  FROM  AMPLITUDE  SPECTRAL  FITTING:  NTS/926/74 
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FIGURE  V-7 

RESULTS  FROM  AMPLITUDE  SPECTRAL  FITTING:  NTS/430/75 
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FIGURE  V-8 

RESULTS  FROM  AMPLITUDE  SPECTRAL  FITTING:  NTS/514/75 
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V-30 


I 

FIGURE  V-8 

j RESULTS  FROM  AMPLITUDE  SPECTRAL  FITTING:  NTS/514/75 

(PAGE  3 OF  3) 
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RESULTS  FROM  AMPLITUDE  SPECTRAL  FITTING:  NTl/603/75 
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FIGURE  V-9 

RESULTS  FROM  AMPLITUDE  SPECTRAL  FITTING:  NTl/603/75 

(PAGE  2 OF  2) 
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(a)  Residual  Distributions 
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the  observed  spectra  and  the  best  fit  theoretical  spectra  is  fairly  good,  espe- 
cially for  the  Rayleigh  waves.  From  the  residual  distributions  shown  in  Fig- 
ures V-6  to  V-21,  the  estimates  of  the  source  parameters  of  each  event  can 
be  represented  in  terms  of  the  probable  range  of  the  estimate  and  the  percen- 
tage frequency  of  that  probable  range.  The  estimations  of  source  parameters 
based  on  the  residual  distributions  are  listed  in  Table  V-8.  To  give  a clearer 
explanation  of  Table  V-8,  consider  the  event  NTl/603/75.  It  has  a depth  esti- 
mate of  0.  5 km  indicated  in  Table  V-7;  while  Table  V-8  says  that  91%  of  the 
members  in  the  sample  set  of  the  tested  source  parameter  (as  described 

in  Subsection  V-A)  indicate  that  the  probable  range  of  the  depth  estimate  is 
from  0.  5 km  to  1 . 5 km. 

To  show  what  the  surface  wave  radiation  patterns  look  like  for 
the  combined  source  of  a point  explosive  source  and  a double-couple  source, 
for  the  event  NTS/1028/5,  Figure  V-22  displays  the  theoretical  Rayleigh  and 
Love  wave  radiation  patterns  from  source  parameter  estimates  given  in  Table 
V-7.  Superimposed  in  this  figure  are  the  observed  radiation  patterns  repre- 
sented by  crosses  (x)  for  the  Rayleigh  wave  and  open  circles  (o)  for  the  Love 
wave.  It  is  hard  to  determine  radiation  patterns  from  the  observed  values, 
since  there  are  only  eight  Rayleigh  waves  observed  at  each  period  and  only 
four  Love  waves.  Nevertheless,  the  theoretical  radiation  patterns,  in  gen- 
eral, explain  the  observed  values  quite  satisfactorily,  in  the  sense  that  the 
maxima  and  the  minima  in  the  theoretical  radiation  patterns  are  oriented  in 
the  right  directions  to  fit  the  large  and  the  small  observed  values,  respec- 
tively. Also  shown  in  this  figure  are  radiation  patterns  of  the  Love-to- Rayleigh 
(LQ/LR)  ratios.  The  theoretical  LQ/LR  radiation  pattern  agrees  quite  well 
with  the  four  observed  values  at  most  of  the  periods.  This  indicates  that  the 
travel  path  attenuation  correction  using  Tryggvason's  attenuation  curve  prob- 
ably is  fairly  reasonable. 
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ESTIXLATIONS  OF  SOURCE  PARAMETERS  OBTAINED  BY  AMPLITUDE 
SPECTRAL  FITTING  BASED  ON  THE  RESIDUAL  DISTRIBUTIONS: 
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The  source  parameter  estimates  of  these  two  events  are  obtained  by  using 
observation  spectra  at  the  SDCS  stations. 
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2. 


The  EKZ  Events 


For  the  selected  EKZ  events,  estimations  of  source  parameters 
based  on  the  minimum- residual  criterion  are  listed  in  Table  V-9.  Those  based 
on  the  residual  distributions  are  listed  in  Table  V-lO.  The  residual  distribu- 
tions with  respect  to  the  estimated  depth,  dip  angle,  slip  angle,  and  strike 
angle  are  shown  in  Figures  V-23  to  V-30,  together  with  the  spectral  fit  for 
each  event.  Except  for  the  events  EKZ/1214/3  and  EKZ/1029/5,  the  source 
parameter  estimates  of  the  selected  EKZ  events  suggest  that  the  depths  of 
these  events  are  all  very  shallow,  0.  5 km.  The  double-couple  sources  which 
model  tectonic  strain  releases  associated  with  these  events  have  a purely  dip- 
slip  type  of  mechanism,  and  the  strike  directions  vary  among  these  events. 

The  F values  vary  from  1, 0 to  1.5  among  these  events.  Judging  from  the 
estimated  F values  of  these  EKZ  events,  the  double-couple  source  is  par- 
tially responsible  for  the  observed  surface  waves  of  these  events,  although 
the  exact  mechanisms  (earthquake  triggering,  directional  cracking,  or  cavity 
distortion)  whicn  release  the  tectonic  strain  energy  associated  with  these 
events  are  not  known. 

The  theoretical  surface  wave  radiation  patterns,  together  with 
the  observed  values  (indicated  by  the  crosses  (x)),  generated  by  using  the 
source  parameter  estimates  listed  in  Table  V-9,  are  shown  in  Figure  V-31 
for  the  event  EKZ/704/76.  This  event  has  the  largest  number  of  observed 
values  among  these  selected  EKZ  events.  Again,  as  for  the  event  NTS/1028/5, 
it  is  hard  to  form  the  observed  radiation  patterns  from  only  10  observed  values 
of  the  Rayleigh  wave  and  4 observed  values  of  the  Love  wave  at  each  period. 
Nevertheless,  the  theoretical  radiation  patterns  fit  most  of  the  observations. 
The  theoretical  LQ/LR  radiation  patterns  explain  the  observed  values  very 
well.  Again,  this  indicates  that  Tryggvason's  attenuation  curve  is  fairly 
representative  for  the  average  surface  wave  attenuation  along  the  travel  paths 
from  EKZ  to  the  available  stations. 
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3. 


The  PNE  Events 


For  the  selected  PNE  events,  estimates  of  source  parameters 
based  on  the  minimum- residual  criterion  are  listed  in  Table  V-11.  Those 
based  on  the  residual  distributions  are  listed  in  Table  V-12.  Again,  as  ob- 
served for  the  NTS  and  EKZ  events,  the  depth  estimates  of  these  events  are 
all  very  shallow  (0.  5 km).  The  dip,  slip,  and  strike  angle  estimates  indicate 
that  the  double-couple  sources  which  model  the  tectonic  strain  releases  asso- 
ciated with  the  events  PNE/820/72  and  PNE/729/76  probably  have  very  similar 
orientations,  i.  e.  , purely  dip- slip  and  striking  close  to  the  E- W direction.  The 
estimated  F values  for  the  three  PNE  events  range  from  1.0  to  1.25,  indicat- 
ing the  existence  of  the  double-couple  component.  The  residual  distributions 
with  respect  to  the  estimated  source  parameters  are  shown  in  Figures  V-32 
to  V-34,  together  with  the  spectral  fit  for  each  event.  The  spectral  fits  for 
these  events  are  thought  to  be  fair. 

For  the  sake  of  completeness,  theoretical  surface  wave  radia- 
tion patterns  generated  by  using  the  source  parameter  estimates  listed  in 
Table  V-11  are  shown  in  Figure  V-35  for  the  event  PNE/729/76.  Again  5 
observed  Rayleigh  wave  values  and  3 observed  Love  wave  values  available 
at  each  period  are  also  plotted  in  this  figure.  The  general  agreement  between 
the  theoretical  Rayleigh  and  Love  wave  radiation  patterns  and  the  observed 
values  is  good.  Especially,  the  three  observed  LQ/LR  values  fall  quite  close 
to  the  theoretical  LQ/LR  radiation  patterns  for  most  of  the  periods. 

4.  Comparison  Among  the  Selected  NTS,  EKZ,  and  PNE  Events 

From  source  parameter  estimates  of  the  selected  NTS,  EKZ 
and  PNE  events  discussed  above,  it  seems  that  differences  in  the  source 
parameter  estimates  are  not  obvious  among  these  three  groups  of  explosion 
events.  The  results  of  comparing  these  three  groups  of  explosion  events 
can  be  interpreted  as  follows: 


V-124 


ESTIMATIONS  OF  SOURCE  PARAMETERS  OBTAINED  BY  AMPLITUDE 
SPECTRAL  FITTING  BASED  ON  THE  RESIDUAL  DISTRIBUTIONS: 


(b)  Spectral  Fits 
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FIGURE  V-33 

RESULTS  FROM  AMPLITUDE  SPECTRAL  FITTING:  PNE/820/72 

(PAGE  2 OF  2) 
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SURFACE  WAVE  RADIATION  PATTERNS  FOR  EVENT  PNE/729/76 
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SURFACE  WAVE  RADIATION  PATTERNS  FOR  EVENT  PNE/729/76 
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SURFACE  WAVE  RADIATION  PATTERNS  FOR  EVENT  PNE/729/76 

(PAGE  8 OF  18) 


PNE/729/76 

DOUBLE  COUPLE  EXPLOSION 
LQ/LR— PERr0D=50*0 


NORTH— 0 DEGREE 


FIGURE  V-35 
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SURFACE  WAVE  RADIATION  PATTERNS  FOR  EVENT  PNE/729/76 
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FIGURE  V-35 

SURFACE  WAVE  RADIATION  PATTERNS  FOR  EVENT  PNE/729/7b 
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• More  good  quality  surface  wave  data  are  available  for  the  NTS 
events  than  for  the  EKZ  or  PNE  events.  There  are  more  Love 
wave  recordings  available  for  the  NTS  events.  However,  lack 
of  good  Love  wave  recordings  for  the  EKZ  events  does  not  nec- 
essarily imply  that  the  EKZ  events  radiated  less  Love  waves. 

In  fact,  for  the  EKZ  events,  if  the  observation  stations  have 
Love  wave  recordings,  the  observed  Love  wave  signals  are 
compatible  with  the  observed  Rayleigh  wave  signals,  in  terms 
of  data  quality  and  signal  level. 

• The  depth  estimates  of  the  selected  NTS,  EKZ,  and  PNE  events 
are  all  very  shallow.  This  is  to  be  expected,  since  these  events 
are  presumed  underground  explosions. 

• The  estimates  of  dip  and  slip  angle  of  the  double-couple  sources 
of  tectonic  strain  release  vary  among  these  events.  Most  NTS 
events  seem  to  have  more  of  a strike-  slip  type  of  mechanism. 
Most  EKZ  events  seem  to  have  more  of  a dip- slip  type  of  mech- 
anism. 

• The  strike  angle  estimates  of  the  double-couple  sources  for 
the  NTS  events  are  generally  in  the  N-S  direction,  while  those 
for  PNE  events  are  close  to  the  E-W  direction.  It  seems  that 
there  is  no  trend  in  the  strike  angle  estimates  of  the  EKZ 
events. 

• The  estimated  F values  vary  among  the  events  in  each  group 
and  among  these  three  groups  of  explosion  events.  They  all 
indicate  that  these  presumed  explosions  have  associated  tec- 
tonic strain  release.  The  range  of  variation  of  the  estimated 
F values  are  not  significantly  different  among  these  three 

groups.  However,  for  the  events  with  the  same  m , the  esti- 

b 

mated  seismic  moments  (given  in  Tables  V-7,  V-9,  and  V-11 


f 

I 
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as  M ) of  the  tectonic  strain  releases  associated  with  the 
selected  NTS  events  are  noticeably  larger  than  those  asso- 
ciated with  selected  EKZ  or  PNE  events.  Also,  the  variations 
of  minimum  residuals,  with  respect  to  the  tested  F values  in 
the  amplitude  spectral  fitting  process,  are  quite  different  be- 
tween the  selected  NTS  events  and  the  selected  EKZ  or  PNE 
events.  Figure  V-36  shows  the  residual  variations  with  respect 
to  the  tested  F values  for  the  events  NTS/1028/5,  EKZ/704/76, 
and  PNE/729/76.  For  the  NTS  events,  the  minimum  residual 
increases  drastically  when  the  tested  F value  increases  from 
the  estimated  F value  and  the  residual  only  increases  very 
slightly  when  the  tested  F value  decreases  from  the  estimated 
F value.  For  the  EKZ  and  PNE  events,  the  residual  increases 
slowly  in  both  directions. 


To  see  if  the  source  parameter  estimates  of  these  selected 

NTS,  EKZ,  and  PNE  events  can  be  used  for  discrimination  among  these 

three  groups  of  events  or  between  earthquakes  and  explosions,  the  estimates 

of  the  source  parameters,  which  are  thought  to  be  relevant,  are  plotted  as 

functions  of  the  event  m^^.  Figure  V-37  displays  the  estimated  F values 

versus  the  event  m,  for  the  selected  NTS,  EKZ,  and  PNE  events.  It  is 
b 

noticed  that  there  is  no  apparent  correlation  between  the  F value  and  the 
event  m^.  That  is,  the  larger  m^  does  not  necessarily  imply  smaller  or 
larger  F values.  The  estimates  of  the  moments  of  the  explosive  sources 
(given  in  Tables  V-7,  V-9,  and  V-11  as  M.^)  for  these  events  as  a function  of 
the  event  m^^  are  shown  in  Figure  V-38.  Also  plotted  in  this  figure  are  the 
previously  obtained  seismic  moment  estimates  of  various  central  Eurasian 
earthquakes  (Tsai  and  Shen,  1972;  Turnbull, et  al.  , 1973,  1974;  Sun,  1976) 
and  the  theoretically  predicted  relation  (given  by  L_  in  the  figure)  between 
the  seismic  moment  and  the  event  for  the  earthquakes  (Tsai,  1972b). 
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The  results  from  the  examination  of  this  figure  can  be  summarized  as  fol- 
lows ; 

• The  theoretically  predicted  relation  between  the  seismic  mo- 
ment and  the  event  m^  for  the  earthquakes  explains  the  pre- 
viously estimated  seismic  moments  of  the  central  Eurasian 
earthquakes  very  well,  for  m^  from  4.  5 to  6.  0. 

• In  general,  for  the  events  with  similar  m^ , the  moment  esti- 

b 

mates  of  the  explosive  sources  for  the  NTS  events  are  larger 
than  those  for  the  EKZ  or  the  PNE  events. 

• The  moment  estimates,  , of  the  explosive  sources  for  the 
selected  NTS,  EKZ,  and  PNE  events  can  be  well  fitted  with  a 
straight  line  (indicated  by  L in  the  figure)  parallel  to  the 
theoretically  predicted  line  L^.  for  the  earthquakes.  The  sep- 
aration of  these  two  lines  is  about  one  m^^  unit.  This  separation 
in  the  moment  versus  m^^  plot  can  be  used  to  discriminate  the 
earthquake  population  against  the  explosion  population. 

• Comparing  with  the  conventional  versus  m^  plot,  the  mo- 
ment versus  m^^  plot  seems  to  offer  better  discrimination  in 
several  respects.  First,  the  separation  between  the  earthquake 
populations  is  less  in  the  moment  versus  m^  plot.  Second, 
earthquake  and  explosion  populations  are  less  scattered  in  the 

moment  versus  m plot.  Third,  in  the  moment  versus  m, 
b b 

plot,  the  explosions  at  NTS  and  EKZ  can  be  considered  as  one 

explosion  population;  while  in  the  versus  m^  plot,  the  NTS 

events  and  the  EKZ  events  behave  very  differently,  with  the 

NTS  events  quite  often  extending  into  the  earthquake  population, 

making  the  M versus  m plot  ineffective  for  discrimination 
s o 

purposes. 


V-157 


SECTION  VI 
CONCLUSIONS 


The  far-field  long-period  surface  waves,  both  Rayleigh  and 
Love,  observed  at  the  long-period  seismic  network  have  been  studied  for  the 
presumed  underground  nuclear  explosions  at  the  United  States  Nevada  Test 
Site  (NTS)  and  Russian  eastern  Kazakh  (EKZ),  and  for  the  peaceful  nuclear 
explosions  (PNE)  in  the  north  Caspian  Sea  region,  A total  of  twenty- seven 
recent  NTS,  EKZ,  and  PNE  events  have  been  analyzed  in  terms  of  their  sur- 
face wave  dispersion  characteristics,  the  Rayleigh  wave  attenuation  along  the 
travel  path  involved,  their  Rayleigh  and  Love  wave  spectra,  and  their  seismic 
source  parameter  estimates. 

The  collective  behavior  of  the  Rayleigh  and  Love  wave  funda- 
mental mode  group  velocities  for  each  travel  path  involved  has  been  obtained. 
For  some  travel  paths,  the  theoretical  surface  wave  dispersions  for  the  nor- 
mal continental  or  oceanic  path  can  explain  the  observed  group  velocities 
quite  well.  However,  the  deviations  of  the  observed  group  velocities  from 
the  theoretical  group  velocity  curves  are  very  noticeable  in  many  cases. 
Therefore,  the  collection  of  the  representative  surface  wave  group  velocity 
curves,  which  have  been  obtained  here,  along  the  travel  paths  from  NTS, 

EKZ,  and  the  north  Caspian  Sea  region  to  AI^PA,  NORSAR,  LASA,  VLPE 
stations,  SRO,  and  the  SDCS  stations  will  be  helpful  and  can  be  used  as  guide- 
lines for  the  future  analysis  of  the  surface  waves  of  the  seismic  events  in 
these  three  areas. 

The  Rayleigh  wave  attenuation  coefficients  (for  periods  from 
10  to  50  seconds  at  5 second  increments)  have  been  obtained  for  the  continen- 
tal and  mixed  paths  from  NTS  to  the  observation  stations,  for  those  from 
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EKZ  to  the  observation  stations,  for  oceanic  paths  from  the  NTS  to  stations 
in  the  Pacific  Ocean,  and  for  one  path  across  the  Atlantic  Ocean  (KON-ZLP). 
The  agreement  between  the  attenuation  values  obtained  here  and  those  calcu- 
lated by  other  investigators  is  fairly  good.  In  general,  Tryggvason's  attenua- 
tion curve  has  yielded  quite  reasonable  attenuation  corrections  and  is  thought 
suitable  for  the  average  travel  paths. 

From  the  examination  of  the  observed  Rayleigh  and  Love  wave 
amplitude  spectra  of  these  NTS,  EKZ,  and  PNE  events,  it  has  been  shown 
that  the  seismic  sources  of  these  explosions  have  various  degrees  of  double- 
couple component  in  addition  to  the  explosive  source.  The  evidence  ot  this 
double-couple  component  in  these  explosion  events  is  indicated  by  the  pres- 
ence of  Love  wave  motions  observed  at  various  recording  stations  and  the 
azimuthal  variations  of  the  observed  surface  wave  amplitude  spectra.  How- 
ever, we  have  not  observed  any  case  where  the  Love  wave  amplitude  is  al- 
ways less  than  the  corresponding  Rayleigh  wave  amplitude  at  a given  station, 
as  found  by  Lambert  et  al.  (1972).  The  observed  LQ/LR  spectral  ratios  are 
found  compatible  with  the  average  fundamental  mode  LQ/LR  spectral  ratios 
for  the  shallow  earthquakes  theoretically  modeled  by  double-couple  sources. 
This  does  not  quite  agree  with  the  observations  made  by  Lambert  et  al.  (1972), 
that  the  LQ/LR  spectral  ratios  for  explosions  are  less  than  1.  0 and  are  much 
less  than  those  for  earthquakes.  Although  it  has  been  found  that  there  is  a 
difference  between  earthquakes  and  explosions  in  the  average  frequency  de- 
pendence of  LQ/LR  ratios,  this  difference  is  not  thought  to  be  sufficiently 
significant  to  be  used  as  a discriminative  feature  between  earthquakes  and 
explosions,  as  implicitly  suggested  by  Lambert  et  al.  (1972). 

The  source  parameter  estimates  of  selected  NTS,  EKZ,  and 
PNE  events  have  been  obtained  by  fitting  the  observed  surface  wave  amplitude 
spectra  with  the  theoretically  calculated  surface  wave  spectra  of  a combined 
source,  which  consists  of  the  explosive  source  and  the  double-couple  source. 
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The  double-couple  source  is  used  to  model  the  possible  tectonic  strain  re- 
lease associated  with  the  explosion.  The  depth  estimates  of  these  events  in- 
dicate that  most  of  them  took  place  at  very  shallow  depths  (0.  5 km).  The 
very  shallow  depth  estimates  are  thought  to  be  very  reasonable,  since  those 
events  are  presumed  underground  nuclear  explosions.  The  dip  and  slip  angle 
estimates  of  the  NTS  events  do  not  quite  agree  with  the  assumption  made  by 
some  authors  (Toksoz,  et  al.  , 1965;  Lambert,  et  al.  , 1972;  Mitchell,  1 975) 
that  tectonic  strain  release  associated  with  the  NTS  events  can  be  modeled  by 
a double-couple  source  with  a vertical  strike- slip  type  of  mechanism  (i.  e.  , 

6 = 90°,  A = 0°).  The  strike  angle  estimates  for  the  selected  NTS  events  are 
generally  in  the  N-S  direction,  while  those  for  the  selected  PNE  events  are 
close  to  the  E-W  direction.  However,  there  is  no  general  trend  observed  in 
the  strike  angle  estimates  for  the  EKZ  events.  No  apparent  correlation  be- 
tween the  estimated  F values  and  the  event  m has  been  found  for  these 

b 

events.  The  F values  estimated  for  the  NTS,  EKZ,  and  PNE  events  are 

1.0-1.75,  1.0-1.  5,  and  1.0-1.  25,  respectively.  However,  for  the  events 

with  the  same  m,  , the  estimated  seismic  moments  of  the  tectonic  strain  re- 
b 

leases  associated  with  the  selected  NTS  events  are  in  general  larger  than 
those  associated  with  the  selected  EKZ  or  PNE  events. 

From  the  source  parameter  estimations  for  those  NTS,  EKZ, 
and  PNE  events,  the  moment  versus  m^  plot  has  been  found  to  be  very  prom- 
ising for  the  purpose  of  discrimination  between  earthquakes  and  explosions. 
Several  advantages  of  using  the  moment  versus  m,^  plot  over  the  conventional 
versus  m^  plot  have  been  found,  such  as  less  scattering  in  both  populations 
and  wider  separation,  in  terms  of  m^,  between  the  two  populations.  However, 
it  must  be  mentioned  here  that  the  earthquakes,  shown  in  Figure  V-38  for  the 
moment  versus  m^  plot,  are  modeled  by  the  double-coupls  source  only  instead 
of  the  combined  source  which  is  used  to  model  the  explosion  events  studied 
here.  It  is  believed  that  the  moment  versus  m^  plot  shown  in  Figure  V-38 
will  not  be  significantly  altered  by  modeling  the  earthquakes  with  the  combined 
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source,  since  the  double-couple  source  alone  has  been  producing  very  satis- 
factory results  for  the  earthquakes.  Nevertheless,  for  a fair  comparison, 
the  earthquake  events  should  be  treated  in  the  same  way  as  the  explosion 
events  by  using  the  combined  source  in  the  source  parameter  estimation  for 
earthquakes. 
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